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ABSTRACT 


The  problem  of  estimating  pressure  drops  with  typical  in-line 
bar,  spray  cooler  configurations  is  discussed.  Such  coolers  form 
an  important  component  of  facilities  for  rocket  testing  at  simulated 
altitudes.  In  the  present  report,  the  complex  nature  of  gas-liquid 
flows  is  first  discussed.  This  is  followed  by  a  presentation  of  re¬ 
sults  of  some  simple  experiments  of  air  flow  over  a  tube  bank  in  a 
model  spray-cooler  section  to  determine  the  influence  of  this  parti¬ 
cular  geometry  on  the  drag  coefficient.  Based  upon  the  information 
in  these  two  sections,  a  semi-empirical  analysis  of  two-phase  flow 
in  coolers  is  then  developed  and  some  results  of  the  method  are  pre¬ 
sented.  Finally,  recommendations  are  made  concerning  an  experimental 
program  which  would  enhance  the  further  development  of  prediction 
techniques  in  cooler  design. 
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NOMENCLATURE 

Cross  sectional  area 
Minimum  frontal  flow  area 

Local  width  of  constant  velocity  jet  core,  b^  -  y ^  (Fig.  6) 
Parameter  in  equation  (11) 

Parameter  in  equation  (11) ,  also  width  of  free  shear  layer 
in  jet  analysis  of  section  II 

Half-width  of  space  between  adjacent  spray  bars 

Coefficient  in  equations  (15)  and  (16) 

Heat  capacity 

Hydraulic  diameter 
Pipe  or  duct  diameter 

Diameter  of  tube  in  spray  bank 

Kinetic  energy  rate 

Parameter  defined  in  equation  (38) 

Fanning  friction  factor 

Mass  flux,  mass  flow  per  unit  area 

Gravitational  constant 

Gravitational  acceleration 
Enthalpy 

Specific  enthalpy 

Latent  heat  of  vaporization 

Reference  length 

Parameter  defined  in  equation  (24) 

Exponent  in  equation  (25) 

Mass  flow  rate 

Parameter  defined  in  equation  (25) 

Exponent  in  equation  (15) 

Fluid  static  pressure 


viii 


AEDC-TR- 68-127 


P 

O 

Q 

q 

R 

R 

R 

g 

r 

Re 

s 

TD 

u 

u  . 

3 

V 

V 
w 
X 
x 
x 

V 


O' 

Y 

c 

'll 

9 

X 

u 

P 

T 

W 


Plenum  stagnation  pressure 

Volumetric  flow  rate  in  pipe 
Volumetric  flow  rate  in  jet  analysis 

Holdup  parameter,  average  value  over  the  pipe  length 
Holdup  parameter,  local  value 
Specific  gas  constant 

Pipe  radius 
Reynolds  1  number 

Longitudinal  center-to-center  spacing  of  tubes  in  bank 
Tube  drag 

Local  longitudinal  velocity  in  tube  bank 
Maximum  jet  velocity  in  tube  bank  (Fig.  11) 

Velocity 

Average  velocity  in  pipe 
Weight  flow  rate 

Parameter  defined  in  equation  (21) 

Longitudinal  distance 

Parameter  defined  in  equation  (39) 

Mole  fraction  in  vapor  phase 

Location  of  core  edge  of  free  shear  layer 

Location  of  free  edge  of  shear  layer 

Void  fraction 

Entrainment  parameter,  equation  (51) 

Pressure  loss  coefficient  in  tube  bank 

Dimensionless  width  defined  in  equation  (32) 

Inclination  angle  from  vertical 

Flowing  liquid  volume  fraction  Q  /  (Q  +  6) 

L  Li  G 

Fluid  viscosity 
Fluid  density 

Shearing  stress  at  pipe  wall 
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$  Parameter  defined  in  equation  (17) 

SUBSCRIPTS 

a,b,c  Arbitrary  subscripts  in  equation  (47) 
cw  Cooling  water 

d  Liquid  droplet 

e  Combined  exhaust  and  steam  flows 

eff  Effective 

ex  Exhaust  gas 

B  Location  of  dividing  streamline  of  free  shear  layer 

G  Gas  phase 

i  Integer  index  for  rows  of  spray  bars:  1,2,3,.... 

L , 1  Liquid  phase 

li  Liquid  injected  as  spray 

If  Liquid  flowing,  liquid  film 

m  Mean  value 

n  Integer  index  for  control  volumes:  1,2,3,... 

N  The  last  control  volume 

NS  No  slip,'  homogeneous  flow  (Dukler's  Case  I) 

o  Value  at  minimum  area  between  two  adjacent  tubes  in  bank 

s  Steam 

T  Total  combined  exhaust  and  droplet  flows 

TP  Two  phase 
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SECTION  I 

GENERAL  REMARKS 


1.1  INTRODUCTION 

In  the  performance  testing  of  rocket  propulsion  systems  under 
simulated  high-altitude  conditions,  large  amounts  of  hot  combustion 
gas  must  be  exhausted  from  the  test  cell  and  diffused  to  atmosphe.ric 
pressure.  Considerable  cooling  of  these  gases  is  necessary  in  order 
to  reduce  the  volumetric  flow  rate  as  well  as  to  protect  containing 
surfaces.  A  common  cooling  technique  involves  the  partial  evapora¬ 
tion  of  water  sprays  which,  while  being  both  effective  and  relative¬ 
ly  simple,  results  in  the  further  requirement  of  separating  the  maxi¬ 
mum  possible  amount  of  the  water  vapor  and  droplets  from  the  cooled 
exhaust  gas.  The  presence  of  water  droplets  is  highly  detrimental  to 
rotary  compressors,  and  more  than  the  minimal  amount  of  water  vapor 
in  the  exhaust  requires  extra  compressor  capacity  to  achieve  speci¬ 
fied  performance  criteria. 

The  use  of  water  sprays  to  cool  hot  exhaust  gases  gives  a  simu- 
taneous  flow  of  both  gas  and  liquid  phases  in  the  exhaust  ducts  of 
rocket  test  cells.  Such  two-phase  gas-liquid  flow  is  a  common  occur¬ 
rence  in  many  other  industrial  processes,  with  boiling  and  condensing 
processes  perhaps  being  the  most  familiar  examples  of  two-phase  flow 
processes.  The  transportation  of  oil,  water,  and  gas  simultaneously 
inside  pipes  is  an  integral  part  of  the  production  and  manufacture  of 
petroleum  and  petrochemical  products.  Other  two-phase  flows  which 
might  be  mentioned  are  the  injection  of  liquid  fuels  into  combustion 
chambers  and  the  use  of  liquid  films  to  cool  the  walls  of  rocket  en¬ 
gines  . 

In  all  the  two-phase  flow  processes  mentioned  above,  the  two 
phases  involved  are  gas  and  the  liquid  phases.  Two-phase  flows  could 
involve  the  flow  of  solid  particles  simultaneously  with  either  a 
liquid  or  a  gas,  but  reference  in  this  report  to  "two-phase"  flow  will 
be  restricted  to  flows  involving  only  gas  and  liquid  phases. 

This  report  is  concerned  with  the  prediction  of  pressure  drop 
across  spray  cooler  sections  contained  within  exhaust  ducts  of  rocket 
test  facilities  used  to  evaluate  rocket  performance  at  simulated  high 
altitudes.  In  the  first  portion  of  this  report,  the  current  state  of 
theoretical  approaches  to  the  general  problem  of  two-phase  flow  is  re¬ 
viewed.  These  analyses  presume  fully  developed  flow  conditions  as  a 
simplifying  assumption,  a  condition  not  fulfilled  in  the  operation  of 
these  spray  coolers.  The  spray  coolers  considered  here  are  only  a  few 
duct  diameters  in  length,  and  the  flows  are  thus  confined  to  the 
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developing  flows  characteristic  of  entrance  regions.  However,  the 
summary  of  current  two-phase  flow  theory  will  provide  background 
knowledge  for  consideration  of  the  more  complex  problems  encounter¬ 
ed  in  these  spray  coolers. 

The  discussion  in  Section  II  includes  experimental  results  to 
improve  estimation  of  the  pressure  loss  coefficient  in  subsequent 
expressions  for  the  pressure  drop  associated  with  tube  drag.  The 
necessity  encountered  was  one  of  specializing  a  previous  analysis 
to  the  particular  geometry  of  the  J-5  spray  cooler  located  at  AEDC . 

In  Section  III  a  mathematical  model  for  the  design  of  spray 
coolers  is  presented.  Test  firing  data  were  then  used  to  assess 
certain  features  of  the  model.  Based  on  these  evaluations,  a  cor¬ 
relation  was  obtained  for  the  "droplet  entrainment  parameter"  con¬ 
tained  in  the  model. 

1.2  PRESENT  STATE  OF  KNOWLEDGE 

Because  of  the  large  number  of  publications  which  have  ap¬ 
peared  concerning  two-phase  flow,  a  comprehensive  review  of  the 
literature  will  not  be  presented  in  this  report.  Instead,  a  brief 
survey  of  the  two-phase  flow  literature  will  be  given,  with  des¬ 
criptions  of  particular  work  included  only  if  the  work  represents 
a  significant  advance  or  is  pertinent  as  background  to  the  analysis 
presented  in  this  report.  The  reader  is  referred  to  Ref.  1-6  for 
more  detailed  reviews  of  two-phase  flow,  of  which  Ref.  1-5  also 
contain  extensive  bibliographies  on  the  subject.  References  7-9 
are  entirely  bibliographies  on  two-phase  flows,  and  Ref.  10  gives 
an  extensive  treatment  on  the  injection  and  combustion  of  fuels. 

The  work  by  Gouse  (Ref.  7)  containing  5253  entries  represents  the 
largest  and  most  recent  of  the  bibliographies,  whereas  Kepple  and 
Tung  have  compiled  in  Ref.  8  an  especially  useful  bibliography  of 
2843  entries,  each  containing  an  abstract  of  the  reference. 

The  science  of  fluid  mechanics  has  evolved  mainly  through  the 
study  of  single-phase  flow.  The  effort  expended  over  a  great  many 
years  has  given  an  understanding  of  single-phase  flow  sufficient 
for  analysis  and  design  work  in  many  situations.  Because  many  pro¬ 
cesses  of  economic  and  technical  importance  involve  two-phase  flow, 
there  is  a  similar  motive  to  develop  an  equally  good  understanding 
of  two-phase  flow.  However  the  complexity  of  two-phase  flow,  as 
contrasted  to  single-phase  flow,  has  caused  the  understanding  of 
such  flow  to  be  much  less  complete  than  of  single-phase  flow. 

To  illustrate  the  complexity  of  two-phase  flow  relative  to 
single-phase  flow,  recall  that  the  prediction  of  friction  pressure 
loss  depends  on  five  variables:  Fluid  density,  fluid  viscosity, 
fluid  velocity ,  pipe  diameter,  and  pipe  wall  roughness.  Further¬ 
more,  the  flow  may  be  characterized  as  laminar  or  turbulent,  de¬ 
pending  simply  on  the  magnitude  of  the  Reynolds  number.  However,  the 
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friction  pressure  loss  in  two-phase  flow  has  been  found  to  de¬ 
pend  on  the  viscosity,  density,  and  velocity  of  each  phase,  pipe 
diameter,  pipe  wall  roughness,  surface  tension  at  the  gas-liquid 
interface,  local  acceleration  of  gravity,  and  the  angle  of  incli¬ 
nation  of  the  flow  direction.  Thus,  there  are  eleven  independent 
parameters  upon  which  friction  pressure  drop  depends  in  two-phase 
flow  as  contrasted  to  five  for  single-phase  flow.  The  number  of 
dimensionless  groups  required  to  relate  the  variables  in  two-phase 
flow  is  nine,  whereas  only  three  are  needed  for  single-phase  flow. 
Such  a  large  number  of  independent  variables  makes  it  a  formidable 
task  to  obtain  a  general  correlation  by  dimensional  analysis  and 
experimental  measurements.  Although  certain  dimensionless  groups 
appear  in  most  correlations,  a  generalized  friction  pressure  drop 
correlation  based  primarily  on  dimensional  analysis  has  not  yet 
appeared. 

Although  the  concept  of  a  critical  Reynolds  number  is  suffi¬ 
cient  in  single-phase  flow  to  define  the  transition  from  laminar 
to  turbulent  flow,  the  case  for  two-phase  flow  is  much  different. 
Two-phase  flow  cannot  be  so  simply  classified  into  these  two  flow 
regimes  which  characterize  single-phase  flow  because  the  liquid 
and  gas  phases  are  distributed  over  the  conduit  cross-section  in 
different  ways  depending  on  the  flow  conditions,  and  each  phase  may 
flow  laminar ly  or  turbulently.  Thus,  a  simple  Reynolds  number  cri¬ 
terion  alone  is  insufficient  to  specify  the  type  of  flow  because  it 
does  not  take  into  account  the  phase  distributions  which  occur  for 
different  conditions  of  two-phase  flow.  However,  despite  the  rela¬ 
tive  complexity  of  two-phase  flow  as  compared  to  single-phase  flow, 
investigations  over  the  last  two  decades  have  achieved  improved 
understanding  of  these  flow  phenomena. 

1.3  FLOW  PATTERNS 

The  distributions  of  phases  or  "flow  patterns",  as  mentioned 
above,  are  an  important  consideration  from  a  fundamental  standpoint 
to  a  meaningful  description  of  two-phase  flow  behavior.  Although 
solutions  based  on  mathematical  models  for  a  particular  flow  pattern 
have  not  been  generally  of  much  greater  accuracy  than  those  which 
ignore  the  flow  pattern,  it  is  usually  acknowledged  that  each  flow 
pattern  must  be  analyzed  separately  to  predict  pressure  drop  accur¬ 
ately  in  two-phase  flow.  The  failure  of  flow  pattern  models  to  pre¬ 
dict  friction  pressure  drop  with  greater  accuracy  can  be  attributed 
in  part  to  the  real  difficulty  in  describing  the  flow  structure  in 
detail.  At  present  such  detailed  mathematical  descriptions  of  the 
flow  structure  are  not  available.  Thus,  the  simplifying  assumptions 
and  approximations  made  to  give  a  tractable  mathematical  model  cause 
the  predicted  friction  pressure  drops  to  deviate  from  the  actual 
values . 
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Another  problem  inherent  in  the  use  of  flow  pattern  models 
is  that  of  predicting  what  flow  pattern  will  exist  for  a  given 
set  of  flow  conditions.  No  quantitative  means  for  defining  flow 
patterns  now  exists,  and  flow  patterns  descriptions  are  based  on 
visual  observation  of  the  phase  distributions  through  transparent 
sections  in  a  flow  conduit.  Because  the  results  of  flow  pattern 
experiments  are  reported  as  verbal  descriptions  of  visual  observa¬ 
tions,  descriptions  and  designations  of  flow  patterns  reported  by 
different  investigators  vary.  Moreover,  the  transition  from  one 
flow  pattern  to  another  as  flow  conditions  are  changed  is  gradual, 
not  abrupt.  A  particular  distribution  of  the  phases  may  be  consid¬ 
ered  as  distinct  by  some  investigators  while  others  might  regard 
this  same  distribution  as  a  transition  zone  between  other  flow  pat¬ 
terns  . 

The  flow  patterns  shown  in  Figure  1  represent  a  frequently  used 
set  of  flow  pattern  classifications  for  horizontal  flow.  Because 
the  boundaries  between  patterns  are  difficult  to  define,  some  of  the 
flow  patterns  shown  might  be  considered  merely  as  transition  patterns. 
For  example,  wavy  flow  might  be  considered  as  the  transition  from 
stratified  to  slug  flow.  Spray  flow  could  be  considered  the  limit¬ 
ing  case  of  annular  flow  (very  thin  liquid  film  on  the  pipe  wall) . 
Alternatively,  more  classifications  than  those  shown  in  Figure  1 
have  been  presented;  Knowles  et  al.  (Ref.  11)  proposed  twelve  flow 
patterns  for  horizontal  adiabatic  flow.  Similarly,  one  set  of  flow 
patterns  which  has  been  delineated  for  adiabatic  vertical  upflow  is 
shown  in  Figure  2.  Note  that  the  stratified  and  wavy  patterns  do 
not  exist  for  vertical  flow. 

All  attempts  to  predict  which  flow  pattern  will  exist  for  a 
given  set  of  flow  conditions  have  been  based  on  empirical  correla¬ 
tions  of  experimental  data,  the  most  widely  used  correlation  being 
the  well-known  "Baker  Flow  Pattern  Chart"  for  horizontal  flow  (see 
Fig.  3  and  Ref.  12) .  Notice  that  Baker  used  the  term  "dispersed" 
flow  for  "spray"  flow  (also  called  "mist"  flow  by  some  investiga¬ 
tors)  .  Baker  also  emphasized  that  the  boundaries  between  the  flow 
patterns  are  not  really  lines,  as  shown  on  the  chart,  because  the 
transition  from  one  flow  pattern  to  another  is  gradual.  As  expect¬ 
ed,  the  sequence  of  flow  patterns  which  can  exist  depends  upon  both 
the  gas  flow  rate  and  the  liquid  flow  rate.  A  number  of  other  flow 
pattern  charts  which  have  been  developed  are  also  described  in 
(Ref.  12) . 

A  characteristic  common  to  all  flow  pattern  correlations  is 
that  each  will  predict  flow  patterns  correctly  for  some  flow  con¬ 
ditions,  but  not  for  others.  This  discrepancy  between  predicted 
and  observed  flow  patterns  may  be  attributed  to  the  facts  that  (1) 
flow  patterns  are  classified  according  to  visual  observations,  the 
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boundaries  between  flow  patterns  not  being  well  defined,  (2)  most 
correlations  do  not  include  all  the  pertinent  independent  vari¬ 
ables,  (3)  the  correlations  are  based  on  data  for  only  a  limited 
range  of  flow  conditions. 

To  refine  flow  pattern  predictions  further,  a  quantitative  and 
measurable  parameter  which  characterizes  the  flow  pattern  must  be 
found.  Some  effort  has  been  directed  toward  such  quantitative  flow 
pattern  characterization  (Ref.  13) ,  but  this  work  is  yet  at  an 
early  stage  and  no  general  correlations  have  emerged. 

1.4  PRESSURE  DROP 


The  momentum  integral  equation  for  steady  upward  flow  within 
a  circular  pipe  (inclined  at  an  angle  0  to  the  vertical)  is 


dP 

dx 


8 


d2  g 

p  yc 


_d 

dx 


pDp/2  2 

pv  rdr 


8gCos9 

2 

D  g 

p  c 


fW2 


p  rdr 


(1) 


The  first  term  on  the  right  hand  side  of  equation  (1)  represents 
the  friction  contribution  to  pressure  drop,  the  second  gives  the 
acceleration  contribution,  and  the  third  expresses  the  gravita¬ 
tional  contribution. 


In  one-dimensional  single-phase  flow  the  friction  factor  con¬ 
cept  is  typically  used  to  represent  the  wall  shear  stress  in  the 
friction  pressure  drop  term.  With  the  introduction  of  the  Fanning 
friction  factor. 


f  = 


T  /  (OV 

w 


/  2go, 


the  friction  pressure-drop  term  becomes 
-2pV2f/ (gcDp) 

For  one-dimensional  single-phase  flow  the  acceleration  term  is 

_  (4w) 2  d(l/n) 

,  ^2.2  dx 

nD  )  g 
p  c 

and  the  gravity  term  is 
-  (pg/gc)Cose 


(2) 


(3) 


(4) 


(5) 
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With  these  simplifications  the  total  pressure  drop  is  ex¬ 
pressible  as 


dP  _  -2p y2f 

dx  D  g 
P  C 


(4w) 

/  2,  2 
(ttD  )  g 
P 
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a (i/p) 

dx 


0  g  Cos 
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(6) 


Note  that  cos  0=0  for  horizontal  flow,  and  also  that  the  ac¬ 
celeration  term  is  zero  for  an  incompressible  fluid.  For  a  com¬ 
pressible  fluid  the  density  can  be  expressed  in  terms  of  the  fluid 
pressure  and  temperature  by  an  equation  of  state.  Thus,  in  single 
phase  flow,  the  above  equation  can  be  used  to  determine  friction 
factors  based  on  experimental  measurements  of  total  pressure  drop. 
Conversely,  once  the  friction  factor  has  been  correlated  to  perti¬ 
nent  flow  parameters  in  a  general  manner,  total  pressure  drop  can 
be  calculated  from  the  same  correlation. 

For  two-phase  flow,  the  simplifications  that  were  possible  for 
single-phase  flow  cannot  be  applied.  In  single-phase  flow  the  ac¬ 
celeration  contribution  to  pressure  drop  results  from  the  expansion 
of  a  compressible  fluid  as  the  pressure  decreases  along  the  con¬ 
duit.  In  two-phase  flow  the  gas  phase  expands  as  the  pressure  is 
reduced,  but  the  phase  distribution  also  changes  and  evaporation 
or  condensation  will  change  the  cross-sectional  areas  occupied  by 
the  two  phases,  thus  causing  an  additional  acceleration  contribu¬ 
tion  . 

If  a  two-phase  flow  is  created  by  the  mixing  of  a  liquid  and  a 
gas,  there  is  an  acceleration  contribution  to  the  pressure  drop  as 
the  phase  distribution  undergoes  changes  until  a  quasi-equilibrium 
distribution  is  reached.  The  movement  of  liquid  within  that  phase 
from  a  region  of  low  velocity  to  a  region  of  higher  velocity,  or 
vice-versa,  causes  an  acceleration  of  that  fluid.  The  term  quasi¬ 
equilibrium  is  used  to  acknowledge  that  phase  distribution  is  al¬ 
ways  changing  in  the  axial  direction  (even  in  the  absence  of  eva¬ 
poration,  condensation,  or  entrance  effects)  because  the  phase  dis¬ 
tribution  is  influenced  by  gas  phase  density,  which  is  dependent 
on  pressure.  Thus,  with  the  continuous  decrease  in  pressure  in  the 
flow  direction  there  is  a  concomitant  and  continuous  change  in  the 
phase  distribution.  For  some  flow  conditions,  such  as  those  at  high 
pressures,  small  liquid  flow  rates,  or  in  large  diameter  pipes,  the 
acceleration  contribution  to  the  pressure  drop  may  be  small,  where¬ 
as  the  acceleration  contribution  has  been  shown  to  be  large  (Ref.  14) 
for  flows  at  low  pressures  or  in  small  diameter  pipes. 
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The  gravitational  contribution  to  pressure  drop  in  two- 
phase  flow  can  be  determined  only  if  the  density  of  the  two- 
phase  mixture  is  a  known  function  of  the  radial  coordinate.  Thus, 
the  gravity  contribution,  like  the  acceleration  contribution,  is 
significantly  dependent  upon  the  phase  distribution  or  flow  pattern. 

The  "holdup"  parameter  (R  ) ,  which  is  frequently  encountered  in 
the  two-phase  flow  literature,  is  also  closely  related  to  phase 
distribution.  Holdup  is  defined  as  the  fraction  of  the  total  cross- 
sectional  area  occupied  by  the  liquid  phase  under  actual  flow  con¬ 
ditions.  It  should  be  noted  that  holdup  is  not  simply  equal  to  the 
ratio  of  liquid  volume  flow  rate  to  total  volume  flow  rate,  even  on 

an  average  basis  (R  ) . 

L 

While  total  pressure  drop  measurements  can  be  made  for  two- 
phase  flows,  the  changes  in  phase  distribution  and  local  phase  ve¬ 
locities  required  for  evaluating  the  acceleration  and  gravitational 
contributions  to  pressure  drop  are  extremely  difficult  to  measure. 

For  this  reason,  very  little  of  the  pressure  drop  data  available 
today  for  two-phase  flow  includes  information  which  would  permit 
even  a  very  crude  approximation  of  the  acceleration  contribution. 
Experimental  data  on  liquid  holdup  permit  only  approximate  values 
for  the  acceleration  and  gravity  terms  to  be  calculated. 

Many  of  the  early  experimenters  underestimated  the  magnitude 
of  the  acceleration  term  and  believed  it  to  be  negligible;  conse¬ 
quently,  the  friction  factor  correlations  developed  were  based  on 
total  pressure  drop  instead  of  friction  pressure  drop.  Unless  the 
correlation  between  the  friction  pressure  drop  and  the  flow  varia¬ 
bles  is  the  same  as  that  between  the  acceleration  pressure  drop  and 
the  flow  variables,  a  correlation  of  friction  factor  to  total  pres¬ 
sure  drop  cannot  accurately  predict  total  pressure  drop  under  flow 
conditions  which  produce  acceleration  pressure  drops  different  from 
those  in  the  experiments  from  which  the  correlation  was  obtained. 

Another  factor  not  considered  in  many  of  the  early  experiments 
is  the  "entrance  length"  required  for  the  phase  distribution  to 
reach  a  qua si-equi librium  state.  Acceleration  effects  are  signifi¬ 
cant  in  the  entrance  region,  and  both  acceleration  pressure  drop  and 
friction  pressure  drop  are  functions  of  the  phase  distribution  imme¬ 
diately  after  mixing.  Because  changes  of  phase  distribution  in  the 
entrance  region  depend  on  how  the  two  phases  are  mixed,  the  method 
of  mixing  becomes  an  additional  independent  variable  if  pressure  drop 
measurements  are  taken  too  close  to  the  mixing  point.  Effects  of  the 
method  of  mixing  have  been  observed  as  much  as  220  pipe  diameters 
downstream  of  the  mixing  point  (Ref.  15) . 
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1.5  PRESSURE  DROP  ANALYSES 

Many  experimental  measurements  of  pressure-drop  have  been  made 
to  date,  and  most  investigators  have  presented  equations  based  on 
their  data  to  predict  pressure  drop.  The  many  empirical  correla¬ 
tions  will  not  be  discussed  here  because,  for  the  most  part,  they 
are  reliable  only  for  flow  conditions  identical  to  those  in  the  ex¬ 
periment  from  which  they  were  obtained.  Instead,  four  semi-empiri¬ 
cal  analyses  will  be  discussed  which  typify  different  approaches  to 
the  problem  of  pressure  drop  prediction. 

Homogeneous  Model 

The  study  of  the  fluid  mechanics  of  two-phase  flow  has  ad¬ 
vanced  primarily  through  the  application  and  extension  of  princi¬ 
ples  developed  in  single-phase  flow.  The  homogeneous  model,  which 
is  considered  in  detail  in  Ref.  2,  represents  the  most  direct  appli¬ 
cation  of  single-phase  flow  principles  to  two-phase  flow  in  that  the 
model  assumes  that  both  phases  are  distributed  uniformly  over  the 
pipe  cross-section  and  that  both  phases  flow  with  the  same  velocity  . 
With  these  assumptions  the  two-phase  mixture  is  treated  as  a  single¬ 
phase  fluid  with  average  properties  for  the  mixture.  Although  var¬ 
ious  definitions  for  the  mixture  properties  have  been  proposed,  the 
similarity  analysis  of  Dukler  et  aJL.  (Ref.  16)  shows  that  the  cor¬ 
rect  mixture  properties  consistent  with  the  assumptions  of  the  homo¬ 
geneous  model  are  given  by: 


P 


PI> 


'G 


U  =  U’-jA  +  UG 


(1-x) 

(7) 

(1-x) 

(8) 

Equation  (6)  can  be  used  to  predict  two-phase  pressure  drop 
based  on  the  homogeneous  model  if  the  density  and  viscosity  terms 
are  taken  to  be  the  mixture  density  and  the  mixture  viscosity.  The 
friction  factor  is  correlated  to  the  mixture  Reynolds  number  just 
as  friction  factors  are  correlated  to  Reynolds  number  in  single¬ 
phase  flow.  The  flow  patterns  shown  in  Figures  1  and  2  suggest 
that  the  homogeneous  model  will  be  more  applicable  to  bubble  and 
spray  flows  since  these  phase  distributions  are  more  nearly  homo¬ 
geneous  that  the  others.  In  general,  however,  the  homogeneous 
models  have  been  unsuccessful  in  predicting  pressure  drops  accurate¬ 
ly  over  a  wide  range  of  flow  conditions. 


Lockhart-Martinelli  Model 

The  most  widely  used  correlation  is  that  of  Lockhart  and 
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Martinelli  (Ref.  17) .  Their  model  accommodates  the  nonhomogeneous 
distribution  of  the  two  phases  over  the  cross-sectional  area  and 
also  different  phase  velocities.  The  basic  postulates  upon  which 
the  model  is  based  are  quoted  as: 

1.  "Static  pressure  drop  for  the  liquid  phase  must  equal  the 
static  pressure  drop  for  the  gaseous  phase  regardless  of 
the  flow  pattern,  as  long  as  an  appreciable  radial  static 
pressure  difference  does  not  exist." 

2.  "The  volume  occupied  by  the  liquid  plus  the  volume  occupied 
by  the  gas  at  any  instant  must  equal  the  total  volume  of 
the  pipe." 

Moreover,  each  phase  is  considered  to  flow  separately  with  its  own 
velocity  and  flow  area,  and  it  is  supposed  that  the  pressure  drop 
in  each  phase  can  be  obtained  from  equations  for  single-phase  flow 
using  the  familiar  Fanning  friction  factor  approach. 


Because  the  two  phases  are  considered  to  flow  independently, 
one  can  write 


and 


W„  = 


plVl 


w_  = 


»gagvg 


(9) 

(10) 


where  the  velocities  V  and  V  are  average  velocities  for  the  two 

L  G 

phases.  The  cross-sectional  flow  areas  are  A  and  A  for  the  liquid 

G 

and  gas  phases,  respectively,  and  these  are  related  to  the  "hydrau¬ 
lic  diameters"  DT  and  D_  by 
L»  G 


al  - 

ag  -  K—r- )  (12> 


The  parameters  a  and  b  depend  on  the  phase  distribution. 

Neglecting  acceleration  and  gravitational  contributions, 
pressure  drop  equation  expressed  in  terms  of  either  phase  is 
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2  a  V2  f 

/  dP\  =  _  PG  G  g 
\  dx/  tp  D  g 

va  C 


(14) 


The  friction  factors  are  assumed  to  be  expressible  in  the  Blasius 
form: 
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where  the  values  of  CT  and  n  or  of  C  and  m  depend  upon  whether  the 

L  G 

flow  is  laminar  or  turbulent.  Combining  equations  (9,11,13,  and  15) 
gives 
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Similarly,  combining  equations  (10,12,14,  and  16)  gives 
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Note  that  (dp/dx)  is  just  the  expression  for  pressure  drop 

for  liquid  flow  alone  in  the  pipe,  and  (dp/dx)  is  the  pressure 

G 

drop  for  single-phase  gas  flow. 

Lockhart  and  Martinelli  assumed  that  the  parameters  §  and  s 

Xj  g 

are  only  functions  of  the  parameter 


(21) 


and  used  X  to  correlate  values  of  $  and  §  obtained  from  experi- 

L  G 

mental  pressure  drop  measurements. 

Lockhart  and  Martinelli  postulated  that  the  flow  of  each  phase 
could  be  either  laminar  or  turbulent;  thus,  they  were  able  to  define 
four  flow  mechanisms: 


1)  turbulent  flow  of  both  phases  (turbulent-turbulent) 

2)  laminar  liquid  flow  and  turbulent  gas  flow  (viscous-turbulent) 

3)  turbulent  liquid  flow  and  laminar  gas  flow  (turbulent-viscous) 

4)  laminar  flow  of  both  phases  (viscous-viscous) 


The  flow  of  a  phase  was  arbitrarily  assumed  to  be  laminar  if  the 

Reynolds  number  based  on  the  pipe  diameter  was  less  than  1000  and 

turbulent  if  the  Reynolds  number  was  greater  than  2000  Values  of 

C  ,  C  ,  m  and  n  for  the  various  flow  types  are  given  in /  Table  1. 

L  G 


Figure  4  shows  the  results  of  the  Lockhart-Martinelli  correlation 
for  the  four  cases  noted  above.  These  curves  can  be  used  to  predict 
pressure  drop  correctly  to  within  about  _+  60  precent  over  a  wide 
range  of  flow  conditions,  and  their  model  has  been  used  primarily 
because  of  its  wide  range  of  applicability. 

Certain  limitations  of  the  Lockhart-Martinelli  analysis  should 
be  mentioned.  The  definition  of  the  four  flow  types  is  arbitrary. 
Furthermore,  the  correlation  does  not  account  separately  for  the  ac¬ 
celeration  contribution  to  pressure  drop.  The  assumption  that  $ 

L 

and  §G  are  functions  of  X  alone  has  been  shown  to  be  incorrect,  as 


has  the  assumption  that  the  pressure  drop  for  the  separate  phases  can 
be  expressed  using  single-phase  flow  equations.  However,  in  a  sub¬ 
sequent  work  based  on  an  annular  flow  pattern  and  assumed  velocity 
profiles,  Levy  (Ref.  18)  determined  an  analytical  relationship  be¬ 
tween  $  $  and  X. 

L,  G 


Similarity  Analysis 

Dukler  e_t  a_l.  (Ref.  16)  applied  the  principle  of  similarity  to 
two-phase  flow  and  found  that  correct  expressions  for  the  friction 
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factor  and  Reynolds  number  are: 
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(25) 


Because  the  quantities  M  and  N  as  well  as  R_  and  R„  are  ob1- 

L  G 

viously  not  known,  the  friction  factor  and  Reynolds  number  expressions 
are  of  little  use  in  these  forms.  However,  Dukler  ert  a_l .  obtained 
simpler  expressions  for  friction  factor  and  Reynolds  number  for  two 
cases : 

Case  I:  ,  Local  phase  velocities  are  equal  (no  slip)  and  the  flow 

is  homogeneous.  This  case  is  identical  to  the  homogeneous  flow 
model  discussed  previously.  With  these  assumptions  equations 
(22)  and  (23)  reduce  to 
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where  pNg  and  nNg  are  given  by  equations  (7)  and  (8)  and  1=  Dp . 

Case  II:  Local  phase  velocities  are  not  equal  but  M  and  N  are 

assumed  to  be  unity.  These  assumptions  give 

(28) 


(29) 

In  this  case  R  is  presumed  available  from  a  correlation  for 
L 

holdup. 

For  Case  I,  the  relationship  between  friction  factor  and  Rey¬ 
nolds  number  was  taken  to  be  the  same  as  that  for  single  phase  flow, 
viz . , 

f  =  0.0014  +  0.125  (Re  )_0‘32  (30) 

ns  ns 

For  Case  II,  the  ratio  f^  /  f  (where  f  is  obtained  from  equation 

tp  o  o 

(30)  using  Re^  from  equation  (29)  for  Re  g)  was  correlated  to  \ 

using  friction  pressure  drops  obtained  from  extensive  experimental 
data .  These  friction  pressure  drops  were  obtained  from  total  pres¬ 
sure  drops  values  of  the  acceleration  and  gravitational  contributions 
estimated  from  calculated  holdup  data.  Based  on  a  statistical  com¬ 
parison  of  Case  I,  Case  II,  and  the  Lockhart-Martinelli  analysis, 
Dukler  concluded  that  the  Case  II  correlation  shows  better  agreement 
with  the  data  in  most  cases. 

The  work  of  Dukler  &t  aJL.  is  significant  not  only  because  it  re¬ 
presents  the  first  formal  application  of  similarity  analysis  to  two- 
phase  flow,  but  because  the  Case  II  correlation  is  based  on  over  400 
experimental  data  points  covering  a  wide  range  of  flow  conditions. 

Levy's  Mixing  Length  Analysis 

None  of  the  analyses  mentioned  previously  attempt  to  describe 
the  local  behavior  of  the  two  phases.  Levy  (Ref.  19)  adapted  the 
well-known  mixing  length  concepts  of  turbulent  single-phase  flow  to 
predict  pressure  drop  in  two-phase  flow.  Levy  assumed  that  the  two- 
phase  flow  can  be  represented  by  a  mixture  whose  density  varies  con¬ 
tinuously  along  the  radial  co-ordinate.  The  mass  transfer  component 
of  the  turbulent  shear  stresses  was  retained  in  the  expression  for  the 
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Reynolds  stresses.  The  fluctuating  velocity  and  density  components 
were  related  to  the  average  velocity  and  density  gradients  using  the 
Van  Driest  modification  of  Prandtl's  mixing  length  theory.  The 
mixing  length  for  a  mass  transfer  was  assumed  equal  to  that  for  mo¬ 
mentum  transfer.  Thus,  the  mixing  length  equations  plus  a  potential 
flow  equation  relating  local  time  average  values  of  the  velocity  and 
density  gave  an  ordinary  differential  equation  which  could  be  solved 
numerically  to  obtain  the  density  distribution  over  the  flow  cross- 
sectional  area.  Friction  factors  were  predicted  in  terms  of  a  mean 
density  and  a  Reynolds  number  based  on  total  flow  rate  and  liquid 
viscosity.  The  mean  fluid  density  was  related  empirically  to  known 
flow  variables  using  experimental  data. 

Levy's  mixing  length  model  shows  acceptable  agreement  with  exper¬ 
imental  data  only  for  a  limited  range  of  the  flow  conditions.  However 
the  mixing  length  model  is  an  important  contribution  to  improved 
understanding  of  two-phase  flow  because  it  treats  such  flow  on  a  local 
basis  using  the  conservation  equations.  The  inadequacy  of  the  pres¬ 
sure  drop  predictions  can  be  attributed  to  the  many  simplifying  as¬ 
sumptions  which  were  made  necessary  by  the  present  lack  of  knowledge 
of  these  local  conditions. 

In  summary,  it  may  be  said  that  although  a  certain  understanding 
of  fundamental  aspects  of  relatively  simple  two-phase  flow  behavior 
does  exist,  current  knowledge  is  insufficient  to  permit  accurate  pre¬ 
diction  of  pressure  drops  associated  with  the  concomitant  problems  of 
handling  and  cooling  large  volumes  of  hot  rocket  exhaust  gases  by  use 
of  water  sprays.  Consequently,  other  approaches  described  in  the  sub¬ 
sequent  sections  of  this  report  were  undertaken. 
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SECTION  II 

PRESSURE  LOSS  COEFFICIENT  FOR  TUBE  DRAG 


2.1  INTRODUCTION 

The  problem  of  specifying  the  pressure  drop  coefficient  to  fix 
the  pressure  loss  associated  with  flow  over  the  spray  bar  array  is 
considered  in  this  section.  There  is  apparently  no  previous  work 
sufficient  to  describe  uniquely  the  situation  encountered  within  the 
spray  cooler  problem.  Although  pressure  losses  caused  by  flow  over 
tube  banks  have  been  long  studied  as  part  of  the  design  of  certain 
types  of  heat  exchangers  and  similar  configurations,  none  of  these 
studies  has  considered  the  spray  bar  problem  pertinent  to  this  inves¬ 
tigation.  The  bulk  of  the  available  data  on  pressure  drop  across 
tube  banks  resulted  from  experimental  studies  prior  to  1940,  with 
another  surge  of  information  being  generated  in  the  early  1940's 
because  of  World  War  II.  Typical  of  the  latter  works  are  those  pub¬ 
lished  as  NACA  Wartime  Reports  (Refs.  20  and  21,  for  example). 

The  approach  to  this  part  of  the  spray  cooler  problem  will  be  to 
modify  the  method  given  by  Abramovich  (Ref.  22)  for  predicting  pres¬ 
sure  loss  with  single-phase  flows  over  tube  banks,  and  subsequently, 
to  incorporate  it  in  the  model  to  be  given  for  the  design  of  spray 
cooler  sections.  Some  theoretical  aspects  of  the  analysis  for  the 
tube  drag  associated  with  single-phase  flows  will  be  considered  first. 
Following  this,  results  of  some  experimental  work  done  as  a  part  of 
this  investigation  to  assess  the  validity  of  theory  and  to  specialize 
results  to  the  spray  cooler  geometry  considered  here  will  be  presented. 
Additional  details  of  this  portion  of  the  investigation  may  be  found 
in  Ref.  23. 

2.2  THEORETICAL  CONSIDERATIONS 

A  theoretical  estimate  of  the  static  pressure  loss  across  a 
tube  bank  can  be  obtained  with  a  method  presented  by  Abramovich 
(Ref.  22).  In  this  flow  model  the  pressure  loss  caused  by  the  tube 
bank  is  considered  to  result  from  the  difference  in  the  kinetic 
energy  of  the  flow  between  a  bar  located  at  station  1  in  Fig.  5  and 
a  position  just  ahead  of  the  next  row  of  bars  (station  1').  The  dif¬ 
ference  in  kinetic  energy  appears  as  a  pressure  loss,  because  the 
flow  must  accelerate  to  pass  through  the  next  row  of  bars.  This  loss 
in  total  kinetic  energy  of  the  jet  is  due  to  mixing  of  the  relatively 
high  velocity  fluid  in  the  core  of  the  jet  with  the  low  speed  flow 
which  circulates  in  the  wake  region  behind  the  bars. 

In  determining  the  kinetic  energy  losses,  the  following  assump¬ 
tions  are  made  about  the  flow:  (1)  there  is  uniform  velocity  at 

stations  1  and  2;  (2)  the  region  between  1  and  1'  is  isobaric; 
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(3)  there  is  no  interaction  between  the  flow  passing  between  ad¬ 
jacent  pairs  of  bars;  (4)  the  flow  is  incompressible. 

As  the  fluid  passes  between  a  pair  of  bars,  viscous  effects 
cause  the  jet-like  flow  to  grow  in  width,  and  a  mixing  zone  or 
free  shear  layer  is  formed  at  each  edge  of  the  jet.  As  the  mixing 
regions  spread,  the  average  velocity  of  the  flow  is  reduced.  The 
core  is  eventually  dissipated  by  viscous  action  and  the  two  shear 
layers  merge  at  some  point  B  (See  Fig.  6) . 

The  effective  "edges"  of  the  shear  layer  (See  Fig.  6)  were 
shown  by  Abramovich  to  be  given  by  y^/b  =  0.416  and  y  /b  =  0.584. 

That  is,  the  inner  portion  of  the  mixing  zone  grows  slightly  slower 
than  the  outer  layer.  The  development,  which  is  omitted  for  brevity, 
is  based  on  approximations  to  the  continuity  and  momentum  equations 
(See  Ref.  22,  Chapter  4,  page  148). 

The  thickness  of  the  free  shear  layer  has  been  found  experimen¬ 
tally  to  increase  linearly  with  the  axial  distance,  x,  (See  Fig.  6) 
and  is  given  by  b  =  0.27x.  The  corresponding  velocity  profile  in 
the  mixing  zone  has  also  been  determined  experimentally  and  may  be 
approximated  by  the  polynomial  given  in  Schlichting  (Ref.  24) : 
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u 
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(31) 


where 


(32) 


and  u  is  the  core  velocity, 
o 

To  determine  the  energy  of  the  jet  at  any  longitudinal  station, 
one  must  first  establish  the  boundaries  of  primary  flow  (i.e.,  the 
flow  which  passes  between  the  bars  at  station  1) .  This  boundary  is 
known  as  the  dividing  streamline  since  it  separates  induced  and  pri¬ 
mary  streams.  The  volumetric  flow  rate  for  the  ha If- jet  at  station 
1  is 
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whereas  at  any  longtitudinal  station  the  same  volumetric  flow  rate  is 
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where  An  is  the  width  of  the  constant  velocity  core  and  y  is  the 
I  B 

location  of  the  dividing  streamline.  But  from  Fig.  6  it  is  evident 

that  A.  =  b  -  y. ,  so  that 
1  o  ' 
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Hence,  the  ratio  of  the  flow  rates  yields 
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But  from  equation  (32) 


(33) 
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so  that,  after  substitution  of  equation  (31)  and  the  foregoing 
relation  into  equation  (33)  one  obtains 
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,  as  was  noted  in  the  beginning  y^/b  =  0.416, 
0.8t3*5  -  0.25T|^  =  0.134 
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which  yields  the  location  of  the  dividing  streamline  as 


=  0.515 


(35) 


Now,  with  the  dividing  streamline  known,  the  energy  loss  due 
to  viscous  mixing  may  be  estimated.  At  station  1  the  initial  kinetic 
energy  of  the  jet  is 


p  u3b 
o  o  o 


(36) 


whereas  at  an  arbitrary  location,  the  kinetic  energy  of  that  flow 
between  the  dividing  streamlines  is 
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(38) 


=  1  -  0.092  b/b 


But  b  =  0.27x  so  that 

e'  =  1  -  0.025x 


(39) 


where  x  =  x/b^.  Between  stations  1  and  1'  with  longitudinal  bar 
spacing  S  (See  Fig.  5),  the  ratio  is 


'  =  1  -  0.25  (1  -  „  1  -  0.025  S/bQ 


(40) 
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because  d/2S  «  1  in  the  present  case.  Thus,  the  change  in  the 
kinetic  energy  ratio  is 

0  —  0  1 

Ae'  =  - - -  =  0.025 

S1 

The  pressure  drop  can 
in  the  following  manner. 

.  ,  2 
m  A  u 

C  .  2 

m  u 

o 

From  the  steady-flow  energy  equation  it  is  noted  that 

.  2  . 

AP/p  ~  Au  /g  ,  and  hence, 
c 


S/b  =  0.025  S  (41) 

o 

be  related  to  the  kinetic  energy  loss 
Let  Ae '  =  £  .  Then 


£  =  gcAP/pu^  (42) 

or 

AP  =  £pu2/g  (43) 

o  c 

where 

r  =  0.025  (S/b  )  (44) 

o 

If  the  pressure  gradient  is  uniform  from  row  to  row,  the  pressure 
drop  for  n  rows  of  tubes  may  be  determined  by  multiplying  the  above 
expression  by  n.  This  result  when  combined  with  equation  (41)  gives 


(AP)n 


0 . 025n 


g 


c 


(45) 


The  pUQ  term  will  be  specialized  to  the  two-phase  flow  analysis  in 

Section  III,  but  the  concern  at  this  point  is  whether  or  not  the 
value  of  the  constant  (0.025)  is  appropriate  for  the  spray  cooler 
geometry.  The  analysis  given  above  has  presumed  a  uniform  array  of 
tubes  of  a  constant  diameter.  This  is  not  the  case  for  the  J-5 
cooler,  which  was  used  in  this  work  as  a  typical  or  example  cooler 
configuration.  Thus,  the  assessment  of  the  validity  of  the  theore¬ 
tical  value  of  0.025  for  the  constant  was  necessary  and  done  by 
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recourse  to  the  experimental  work  described  below. 

2.3  EXPERIMENTAL  INVESTIGATION 

A  series  of  experiments  was  made  on  a  scale  model  of  the  AEDC 
J-5  spray  bar  configuration.  The  model  was  constructed  with  wall 
static  taps  at  several  points  along  the  tube  bank  to  determine  the 
effect  of  upstream  constr ictions  on  the  local  pressure  loss  as  well 
as  to  determine  the  overall  pressure  loss.  Also,  provision  was  made 
for  taking  velocity  profiles  both  ahead  of  and  behind  various  rows 
of  bars  to  obtain  a  description  of  the  flow  field. 

Test  data  were  obtained  from  a  tubular  test  cell,  fitted  be¬ 
tween  conical  diverging  and  converging  sections  (See  Fig.  7)  .  The 
test  cell  was  fabricated  from  plexiglass  and  was  secured  between 
the  two  variable-area  sections  by  three  axial  bolts  which  ran  the 
length  of  the  test  section.  Aluminum  plates  on  either  end  of  the 
test  section,  along  with  the  bolts,  held  the  test  section  in  place 
and  served  to  seal  the  chamber  against  leakage.  The  rectangular 
array  of  bars  inside  the  test  section  was  fabricated  from  brass  rod 
and  held  in  place  by  clear  epoxy  cement. 

Static  pressure  taps  were  installed  at  various  stations  along 
the  walls  of  the  test  section.  The  taps  were  made  of  0.050-  inch 
diameter  steel  hypodermic  tubing,  and  were  pressed  into  the  wall 
with  an  interference  fit.  At  the  inside  wall  of  the  test  section, 
the  tap  diameter  was  0.025  inches. 

At  the  same  longitudinal  location  of  each  static  pressure  tap, 
but  displaced  from  it  by  90  ,  were  two  additional  sets  of  holes 
through  the  wall  which  were  provided  to  facilitate  pitot  tube  tra¬ 
verses  directly  behind  a  row  of  bars,  and  directly  in  front  of  the 
adjacent  row  of  bars.  The  pitot  tube  used  for  these  velocity  pro¬ 
file  measurements  was  also  fabricated  from  0.050-inch  diameter 
stainless  steel  hypodermic  tubing.  One  end  of  the  tubing  was  closed, 
and  a  0.025  inch  hole  was  drilled  into  and  perpendicular  to  the 
wall  near  the  midpoint  of  the  tubing. 

The  pitot  tube  was  calibrated  against  a  standard  Prandtl  impact 
tube  over  the  anticipated  range  of  velocities.  As  would  be  ex¬ 
pected,  there  was  some  error  due  to  local  crossflow  effects  when 
the  tube  was  in  a  steep  transverse  velocity  gradient.  However,  this 
was  not  considered  to  be  of  major  consequence  even  though  such  grad¬ 
ients  exist  in  the  flow  field.  Traverses  were  operated  with  a 
micrometer  lead  screw  mechanism  having  forty  threads  per  inch. 

The  test  apparatus  consisted  of  a  plenum  chamber,  test  section, 
pitot  tubes,  and  manometers.  The  test  section  was  preceded  with  a 
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six-foot  length  of  1.25  inch  diameter  copper  pipe  to  assure 
fully  developed,  turbulent,  parallel  flow.  The  plenum  chamber 
was  fed  from  a  receiving  tank  which,  in  turn  was  supplied  by  a 
two-stage  reciprocating  compressor.  The  air  was  passed  through  two 
porous  stone  filters  for  dehumidification  before  it  reached  the 
plenum  chamber.  Maximum  stagnation  pressure  and  temperature  obtain¬ 
able  in  the  chamber  were  100  psig  and  130  F,  respectively.  The  max¬ 
imum  flow  rate  of  this  compressor  is  approximately  1.0  lbm/sec. 

All  pressures  were  recorded  on  two  100-inch  U-tube  mercury 
manometers  which  were  coupled  in  series.  One  leg  of  each  mano¬ 
meter  was  pressurized  by  a  house  compressor  and  the  unknown  pressure 
(probe  or  static)  was  fed  into  the  one  leg,  while  the  fourth  leg 
was  left  open  to  the  atmosphere.  This  technique  permitted  measure¬ 
ment  of  gage  pressures  up  to  200  inches  Hg. 

Flow  rates  were  determined  from  integration  of  measured  veloc¬ 
ity  profiles  which  were  obtained  immediately  upstream  of  the  first 
row  of  bars.  These  flow  rates,  which  covered  the  range  of  compressor 
capacity,  were  plotted  against  the  plenum  chamber  pressure  to  form 
a  calibration  curve  for  the  test  apparatus. 

As  noted  previously,  transverse  profiles  of  velocity  were  ob¬ 
tained  just  after  the  flow  passed  between  a  given  row  of  bars  (i.e., 
immediately  upstream  of  a  static  pressure  tap)  and  just  before  en¬ 
tering  the  reduced  flow  area  at  the  next  row  of  bars.  These  pro¬ 
files  were  taken  approximately  0.025  inches  behind  the  trailing  edge 
of  the  upstream  row  of  bars  and  0.075  inches  ahead  of  the  downstream 
row  of  bars.  Traverses  were  started  at  the  wall  and  readings  were 
taken  at  the  maximum  and  minimum  velocity  points  over  approximately 
the  first  0.60  inches.  Readings  were  taken  at  0.025  inch  increments 
across  the  center  portion  of  the  test  section,  to  a  point  approxi¬ 
mately  0.60  inches  from  the  opposite  wall. 

Because  visual  distortion  through  the  curved  surface  of  the 
plexiglass  tube  prevented  accurate  positioning  of  the  pitot  tube 
at  the  wall,  the  center  spray  bar  was  used  as  a  reference  point  from 
which  transverse  distances  were  measured.  From  the  center  bar  the 
pitot  tube  was  translated  by  an  amount  equal  to  the  test  section  ra¬ 
dius  and  this  position  was  taken  as  the  effective  wall  location.  The 
inherent  inaccuracy  in  this  method  and  the  slack  in  the  lead  screw 
mechanism  made  a  slight  translation  of  the  velocity  distributions 
necessary  when  the  profiles  were  plotted. 

2.4  DISCUSSION  OF  RESULTS 

Typical  velocity  profiles  found  at  various  stations  in  the  test 
cell  are  shown  in  Fig.  8.  Also  indicated  on  this  plot  is  the  rela¬ 
tive  position  of  the  bars  in  the  flow  field.  The  base  line  of  the 
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velocity  profile  corresponds  to  the  location  of  the  probe.  It  is 
observed  that  the  profiles  exhibit  the  expected  jet-like  peaks  just 
downstream  from  each  row  of  bars.  The  actual  distributions  are 
seen  to  be  nearly  parabolic  initially  and  exhibit  considerable  decay 
(or  flattening)  before  reaching  the  next  row  of  bars.  A  very  slight 
decrease  in  the  successive  values  of  maximum  velocity  is  also  evident. 

The  velocity  profile  upstream  from  each  row  exhibits  the  usual 
shape  of  a  slightly  retarded  (decelerated)  pipe  flow.  It  also  indi¬ 
cates  that  no  flow  separation  occured  in  the  diverging  portions  of 
the  test  cell  upstream  of  the  first  row  of  tubes.  Crossflow  effects 
on  the  pitot  tube  probably  cause  the  measured  velocities  in  the  steep 
gradient  regions  near  the  edges  of  the  jets  to  be  somewhat  higher 
than  the  actual  speeds.  On  the  other  hand,  velocities  near  the  jet 
peaks  should  be  much  more  accurate. 

Wall  pressure  distributions  for  four  mass  flow  rates  (See  Table 
I)  are  shown  in  Fig.  9.  The  wall  pressures  have  been  normalized  with 
the  plenum  (or  receiver)  pressure,  P  .  This  latter  pressure,  it 
should  be  noted,  is  considerably  larger  than  the  local  stagnation 
pressure. 


Values  of  P  may  be  obtained  from  Table  I  for  the  four  runs. 

The  number  of  rows  of  bars  located  between  successive  taps  is  also 
indicated  in  Fig.  9.  It  is  evident  that  the  wall  pressure  distribu¬ 
tions  were  very  nearly  linear,  as  was  postulated  in  the  foregoing 
theoretical  model.  It  should  be  noted  that  the  last  static  tap  was 
located  between  the  seventh  and  eighth  rows  of  bars  so  that  the 
gas  expansion  loss  which  occurs  downstream  of  the  last  row  of  tubes 
would  be  excluded,  thus  affording  a  better  comparison  with  theory. 

The  dimensionless  pressure  gradient  shown  as  a  function  of  the 
mass  flux  in  Fig.  10  increases  with  the  flow  rate,  as  would  be  ex¬ 
pected.  This  information  may  be  used  to  estimate  pressure  drops  for 
other  flow  conditions  with  similar  geometries. 

In  determining  experimental  values  of  the  pressure  drop  coeffi¬ 
cient,  £ ,  it  was  found  that  equation  (42)  was  somewhat  undesirable 
because  of  its  dependence  upon  density,  this  being  difficult  to 
infer  from  other  measured  quantities.  However,  a  more  tractable 
expression  for  £  in  terms  of  easily  accessible  quantities  is 


C 


=  g^AP/np 


2 

u 
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=  g  A P/nGu 
c  o 


(46) 


where  uq  is  the  average  velocity  at  the  minimum  area,  and  G  is 
mass  flux,  also  based  on  the  minimum  area. 
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Because  velocity  profiles  were  obtained  just  downstream  of  a 
row  of  bars,  i.e.,  displaced  downstream  from  the  minimum  area  by  a 
distanceequal  to  one-half  the  bar  diameter  plus  the  clearance  be¬ 
tween  the  bar  and  the  pitot  tube  (a  total  of  0.070  in.),  it  was 
necessary  to  calculate  the  average  velocity,  U  ,  from  the  measured 
profiles.  This  calculation  was  based  on  a  parabolic  velocity  dis¬ 
tribution  at  the  probe.  The  relation  between  and  the  maximum 


j  et 


velocity 


U  . 
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(See  Fig.  11)  is 


2b  u  =  (2/3)  U,  (2b  +d) 

o  o  3  o 

or,  for  the  J-5  cooler  geometry, 


u  =  1.09  u, 
o  3 


The  ratio  of  theoretical  to  experimental  values  of  £  is  shown 
in  Fig.  12.  The  upper  curve  incorporates  the  theoretical  predic¬ 
tion  of  equation  (46)  which  is  valid  for  s/bQ  =  7.23.  It  is  ob¬ 
served  that  the  theory  tends  to  overestimate  the  pressure  loss  for 
three  out  of  the  four  runs.  Results  of  an  analysis  similar  to  that 
in  Section  2.2,  but  including  variable  density,  demonstrated  that 
compressibility  effects  were  indeed  negligible,  as  has  been  assumed. 

Further  consideration  of  possible  sources  of  error  lead  to  an 
exploration  of  the  effects  of  multiple  bar  diameters  in  the  test 
model  (See  Fig.  8) .  The  smaller-diameter  bars  near  the  wall  result 
in  a  change  in  s/bQ  from  7.23  to  6.23  for  that  region  between  the 

two  sizes  of  bars.  Furthermore,  for  the  region  between  the  smaller 
bars  and  the  wall,  the  equivalent  s/b  is  3.15.  The  predicted 

value  of  £ ,  based  on  these  smaller  s/b^  values,  decreased  from  an 

intitial  0.180  to  0.156  and  0.079,  respectively.  An  effective  value 
of  £  for  the  complete  test  section  was  approximated  by  prorating 
the  three  £  values  on  a  cross  sectional  area  basis*.  This  yields 
the  relation 


(47) 


=  (0.180)  (0.584)  +  (0.156)  (0.270)  +  (0 . 079)  (0 . 145) =0 . 158 
The  effective  £  was  thus  14  per  cent  smaller  than  that  given  by 

*The  wall  shear  stress  was  neglected  after  preliminary  estimates 
confirmed  that  its  contribution  was  indeed  negligible. 
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equation  (41) .  The  resulting  ratio  of  theoretical  to  experimental 
values  of  £  is  shown  as  the  lower  curve  of  Fig.  12.  It  is  seen 
that  the  maximum  error  between  theory  and  experiment  has  been  re¬ 
duced  from  33  to  18  percent. 

Figure  13  shows  a  comparison  of  the  present  results  with  those 
of  Abramovich  (Ref.  22)  and  Pierson  (Ref.  25)  for  in-line  arrays  of 
bars.  The  theoretical  pressure  loss  coefficient  is  also  shown.  More¬ 
over,  in  this  figure  are  the  G  values  for  the  present  study  and  for 
Pierson's  data.  The  corresponding  information  was  not  known  for  the 
Abramovich  study.  It  is  noted  that  the  data  of  Pierson  was  obtained 
at  a  flow  rate  and  pipe  Reynolds  number  considerably  lower  than  in 
the  present  study.  This  may  account  for  a  portion  of  the  disagree¬ 
ment  between  the  theory  and  measured  values. 

Based  on  these  results,  the  theoretical  value  of  0.025  for  the 
constant  in  the  pressure  loss  coefficient  is  considered  to  be  too 
large  for  the  J-5  spray  cooler  geometry.  Consequently,  consistent 
with  the  concept  given  in  equation  (47)  and  the  behavior  noted  in 
Fig.  12,  the  value  of  the  constant  will  be  taken  as  approximately 
0.020.  Thus  the  pressure  loss  coefficient  is  taken  as 

C  =  0.020  (S/b  )  (48) 

o 

for  the  J-5  geometry,  with  b^  being  the  half  spacing  between  the 

larger  diameter  spray  bars.  It  should  be  emphasized  that  equation 
(48)  is  unique  for  the  J-5  configuration  (See  Fig.  14);  however, 
either  of  equations  (48)  or  (44)  will  give  a  reasonable  estimate 
of  the  pressure  loss  coefficient  for  spray  bars  arrayed  in  a  typical 
in-line  configuration. 
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SECTION  III 

APPROXIMATE  ANALYSIS  OF  SPRAY  COOLERS 


3.1  THE  PROBLEM 

The  design  problem  can  be  stated  as  follows:  given  exhaust 
flow  parameters  at  the  inlet  of  the  spray  cooler  section,  cooling 
water  parameters,  and  specified  spray-bar  configuration  and  duct 
geometry,  determine  the  number  of  rows  of  spray  bars  required  to 
cool  the  exhaust  stream  to  a  temperature  level  suitable  for  the 
subsequent  compression  process  and  predict  the  pressure  drop  which 
is  incurred  by  flow  through  the  spray  bank.  Moreover,  it  is  de¬ 
sirable  to  minimize,  if  possible,  the  amount  of  water  droplet 
carryover  to  the  compressors.  With  a  predictive  model  for  spray- 
cooler  design,  parametric  evaluations  are  then  possible  to  suggest 
optimal  operating  conditions  and  design  features. 

Although  large  amounts  of  cooling  water  are  sprayed  into  the 
cooler,  only  a  fraction  of  the  cooling  water  will  augment  the  ex¬ 
haust  gas  flow  as  steam  and  droplets.  Indeed,  with  large  values 
for  the  ratio  of  mass  flow  rate  of  injected  water  spray  to  exhaust 
gas  flow  rate  (typically  about  50:1),  most  of  the  injected  water 
will  fall  out  and  flow  along  the  bottom  of  the  duct  to  drain  sumps 
beyond  the  cooler  section.  Thus,  the  cross-sectional  area  of  the 
duct  occupied  by  the  two  phase  flow  will  not  be  circular  since  the 
liquid  will  form  only  a  thin  film  on  the  upper  part  of  the  duct 
surface  and  most  of  the  water  will  fill  the  lower  portion  of  the 
duct.  This  behavior  is  only  qualitatively  described  because  it  will 
depend  upon  mass  flow  rates,  duct  geometry,  and  the  configuration 
of  drain  sumps. 

Thus,  the  phase  distribution  within  the  cooler  section  pro¬ 
bably  will  be  a  combined  form  of  annular  mist  flow  and  stratified 
wave  flow,  the  dispensed  droplet  -  exhaust  gas  flow  being  in  the 
annular  mist  flow  regime  and  the  excess  cooling  water  flow  along 
the  bottom  of  the  duct  being  in  the  wave  flow  regime.  Moreover, 
as  mentioned  in  Section  1.1,  the  flow  within  the  cooler  section  is 
not  fully  developed  in  any  sense  because  of  the  small  length-to- 
diameter  ratio  of  the  spray  cooler  and  the  complex  nature  of  the 
two-phase  flow,  especially  since  the  water  is  injected  at  each 
successive  row  of  spray  bars.  Thermodynamic  equilibrium  between 
the  phases  flowing  under  these  conditions  probably  does  not  exist, 
although  it  is  reasonable  to  suppose  the  exhaust  gas  flow  leaving 
the  cooler  section  will  be  saturated  with  steam  because  of  the  large 
ratios  of  water-spray  to  exhaust-gas  flow  rates  used.  Furthermore, 
the  entrained  droplet  flow  is  probably  not  in  mechanical  equilibrium 
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with  the  gas  phase  flow,  the  droplets  being  accelerated  in  typical 
operation  of  the  J-5  spray  cooler. 

These  features  of  the  problem  and  the  current  state  of  know¬ 
ledge  about  the  simultaneous  phenomena  involved  augur  against  a 
precise  description  or  model  for  spray  cooler  performance.  Thus, 
the  description  given  here  must  necessarily  be  an  oversimplification 
and  it  represents  only  a  first  attempt  to  obtain  a  model  suitable 
for  the  design  of  spray  coolers. 

3.2  THE  DESIGN  MODEL 

The  model  to  be  presented  is  comprised  of  solutions  to  the 
momentum  and  energy  equations  augmented  by  arbitrary  assumptions 
and  expressions  for  mass  conservations.  The  energy  equation  is 
used  to  calculate  the  temperature  of  the  exhaust  flow  leaving  the 
cooler  section,  this  temperature  being  the  major  process  desigr 
criterion.  The  momentum  equation  is  used  to  obtain  the  outlet 
pressure,  with  which  the  pressure  drop  through  the  cooler  is  cal¬ 
culable.  The  solution  to  the  problem  is  complicated  by  the  coupled 
nature  of  the  momentum  and  energy  equations  and  by  the  discrete 
conditions  caused  by  the  arrangement  of  successive  rows  of  spray 
bars  in  the  tube  bank. 

Although  it  is  possible  to  write  the  energy,  mass,  and  momen¬ 
tum  conservation  equations  to  describe  the  phenomena  involved  in 
the  spray  cooling  process,  the  present  state  of  knowledge  for  two 
phase  flow  requires  that  a  number  of  simplifying  assumptions  be 
made  if  useful  expressions  are  to  be  derived  from  the  basic  equa¬ 
tions.  The  initial  assumptions  made  in  the  model  given  here  in¬ 
clude  the  following:  (1)  ideal  gas  behavior  is  assumed  for  the 
steam  and  exhaust  gas  components  (2)  frictional  losses  at  the  duct 
wall  are  neglected  (3)  momentum  losses  incurred  in  the  acceleration 
of  the  droplets  to  the  exhaust  flow  velocity  are  neglected  (4) 
momentum  exchanges  as  droplets  enter  the  liquid  film  or  are  entrained 
from  the  liquid  surface  are  ignored. 

The  first  assumption  given  above  is  probably  not  too  severe 
because  pressures  are  reasonably  low  and  temperatures  are  high  in 
much  of  the  cooler  section.  The  restrictions  inherent  in  the  as¬ 
sumed  ideal  behavior  may,  moreover,  be  alleviated  by  use  of  a  better 
equation  of  state  for  the  exhaust  mixture.  However,  the  added  com¬ 
plexity  seems  unwarranted  at  this  stage  of  the  theoretical  develop¬ 
ment.  The  neglect  of  frictional  losses  at  the  duct  wall,  taken  as 
the  second  assumptions,  is  likewise  probably  reasonable.  The  major 
momentum  effects  will  be  those  associated  with  tube  drag  and  with 
the  injected  water  spray.  Furthermore,  the  frictional  losses  at  the 
wall  are  difficult  to  estimate  under  the  best  of  two-phase  flow  con¬ 
ditions,  as  was  discussed  in  Section  I;  the  flow  conditions  within 
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the  spray  cooler  make  such  estimation  even  more  formidable. 

The  third  and  fourth  assumptions  are  less  tenable  than  the 
former  assumptions.  There  is  obviously  momentum  exchange  in  terms 
of  droplet  drag  as  these  particles  are  accelerated  to  the  exhaust 
flow  velocity.  Furthermore,  the  problem  is  additionally  complicated 
with  mass  transfer  by  condensation  or  evaporation  of  droplets.  The 
momentum  exchanges  between  the  droplets  and  liquid  film  is  even  a 
more  complex  problem.  The  interfacial  conditions  at  the  liquid 
film  surface  are  not  amenable  to  any  simple  description.  Much  more 
detailed  knowledge  of  these  features  of  the  flow  must  be  gained 
before  they  can  be  treated  realistically  in  a  model,  such  as  that 
presented  here. 

The  extent  to  which  the  model  is  deficient  for  taking  such  as¬ 
sumptions  depends  upon  the  relative  magnitudes  of  these  effects 
which  are  ignored  and  those  which  are  retained  in  the  model.  It  is 
believed  for  the  spray  cooler  problem  that  these  neglected  effects, 
although  important,  are  not  of  major  importance  when  compared  to 
other  momentum  exchanges  which  are  considered. 

The  analysis  is  based  on  the  configuration  and  control  volumes 
shown  in  Fig.  15,  different  control  volumes  being  selected  for  the 
estimation  of  losses  caused  by  tube  drag  and  by  other  mechanisms. 

It  is  assumed  that  the  inlet  pressure  P  ,  inlet  temperature  T  , 
mass  flow  rates  of  rocket  exhaust  and  s?eam  injected  upstream°from 
the  spray  cooler,  and  the  rocket  exhaust  gas  composition  are  known. 
The  steam  formed  as  a  combustion  product  in  the  exhaust  is  added  to 
the  steam  flow  rate  from  the  ejector  to  obtain  the  total  steam  flow 
rate  entering  the  spray  cooler.  Moreover,  any  condensable  components 
of  the  rocket  exhaust  (such  as  alumina  in  the  case  of  solid  propel¬ 
lants)  are  also  subtracted  from  the  exhaust  gas  flow,  and  the  exhaust 
gas  properties  and  flow  rates  are  then  calculated  on  the  basis  of  the 
"noncondensable"  exhaust  flow.  Thus,  this  "exhaust  gas"  flow  rate 
W  is  constant  through  the  spray  cooler. 

The  cooling  water  W  ^  sprayed  into  a  control  volume  may:  1)  be 

entrained  as  droplet  flow  with  the  combined  exhaust  and  steam  flow 
W^,  2)  be  vaporized  to  augment  the  steam  flow  W  ,  or  3)  fall  out  to 

increase  the  liquid  flow  W  .  The  3roPl-et  flow  entering  the  control 

volume  may  leave  the  volume  as  droplet  flow,  be  vaporized  and  leave 
as  steam,  or  fall  into  the  liquid  flow.  The  steam  flow  into  the  con¬ 
trol  volume  will  be  assumed  either  to  leave  as  steam  or  to  condense 
out  to  increase  the  liquid  flow  (that  is,  condensation  is  assumed  to 
occur  only  on  the  peripheral  liquid  film  surface  rather  than  on  the 
droplets  or  to  form  additional  droplets.)  The  liquid  flow  entering 
the  control  volume  is  assumed  to  leave  the  control  volume  in  the  same 
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manner;  that  is,  the  net  transfer  of  water  droplets  between  the 
liquid  and  the  core  flow  is  to  the  liquid  surface.  Moreover,  drain 
sumps  are  assumed  not  to  be  located  within  the  spray  cooler  section. 

Mass  balances  written  for  the  n-th  control  volume  give 

W  ,  =  W  =  W  (49) 

exn-1  exn  ex 


for  the  exhaust  gas  flow  (as  defined  above)  and 


W,  _  ,  +  W  .  +W,  ,  +W.  .  =  W.  ,  +W  +W . 

lfn-1  sn-1  dn-1  li  lfn  sn  dn 


(50) 


for  the  water. 


At  this  point  in  the  analysis  it  is  necessary  to  postulate 
that  the  cooling  spray  flow  rate  is  sufficiently  large  that  the 
entrained  droplet  flux  will  depend  only  upon  the  "carrying  capacity" 
of  the  total  exhaust  stream.  It  is  assumed  for  this  purpose  that 
the  droplet  flow  rate  is  proportional  to  the  total  exhaust  flow  and 
may  be  expressed  as 


dn 


=  y  (g  + 
n  sn 


G  ) 
exn 


or 


W„  =  y  (W 
dn  Tn  sn 


+  W  ) 
exn 


(51) 


which  is  a  definition  of  y  ,  an  "entrainment  parameter".  While 

n 

such  a  definition  of  a  droplet  entrainment  parameter  might  be  con¬ 
sidered  prestidigitation,  and  while  such  a  definition  admittedly 
ignores  all  description  of  droplet  entrainment  mechanisms,  the  pre¬ 
sumption  of  such  a  postulate  may  be  almost  ineluctable  on  a  simple 
phenomenological  basis,  consistent  with  the  assumption  stated  above. 


Thus,  the  change  in  the  liquid  flow  rate  across  the  control  vol¬ 
ume  can  be  expressed  in  terms  of  the  injected  spray  flow  rate,  the 
steam  flow  rates,  and  the  exhaust  gas  flow  rate  using  equations  (49) , 
(50) ,  and  (51)  to  give 


w.  _  -  w.  =  w_ ,-w 

lfn  lfn-1  li  sn 


+W  -y  (W  +W  )4y  ,  (W  .  +W  )  (52) 

sn-1  n  sn  ex  'n-1  sn-1  ex 


It  will  be  convenient  in  the  development  of  the  engineering 
model  for  the  prediction  of  pressure  drop  in  spray  coolers  to  con¬ 
sider  separately  that  pressure  loss  caused  by  tube  drag  and  that 
associated  with  other  momentum  losses.  For  typical  operating  condi¬ 
tion,  these  two  effects  will  normally  oppose  each  other,  the  tube 
drag  causing  a  pressure  loss  and  the  injected  water  spray  being 
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sufficient  to  cause  a  pressure  rise  through  the  cooler  section. 

The  momentum  losses  other  than  those  caused  by  tube  drag  are 
expressed  by  a  momentum  balance  written  across  the  n-th  control 
volume  as 


g  AAP 
c 


=  (w  +w  +wj  V  -  {w  +w  +wj  _V  , 
ex  s  d  n  n  ex  s  d  n-1  n-1 


-W  .V  +W  V  -  W  V 

li  li  If n  If n  lfn-1  lfn-1 


(53) 


in  which  the  simplifying  assumptions  mentioned  earlier  pertaining 
to  momentum  losses  have  been  taken.  If  it  is  assumed  that  the 
liquid  flows  along  the  duct  wall  at  a  constant  velocity  V  ,  and 

using  equation  (52) ,  one  can  express  equation  (53)  as 


g  AAP 1  =  (1+v  )  (W  +W  )  V  -  ( 1+v  J(W  +W  JV  -W,  .V.  . 

c  n  ex  sn  n  n-1  ex  sn-1  n-1  li  li 

+W  -v  (W  +W  )+y  (W  +W  _)]v.,-  (54) 

L  li  sn  sn-1  ^n  ex  sn  T  n  ex  sn-1  J  If 


n-1 


The  pressure  loss  caused  by  tube  bank  drag  is  predicted  using 
an  analysis  patterned  after  that  by  Abramovich  (Ref.  22)  and  re¬ 
sults  of  the  work  described  in  Section  II.  However,  this  analysis 
incorporates  modifications  to  account  for  the  two  phase  flow, 
whereas  the  Abramovich  analysis  applies  only  to  single-phase  flow 
over  an  array  of  tubes. 


The  modifications  are  based  on  the  concept  of  the  "homogeneous 
model"  (see  Section  1.5)  in  which  it  is  assumed  that  the  liquid 
droplets  and  exhaust  gases  flow  as  a  homogeneous  mixture  with  zero 
relative  velocity  between  the  phases.  In  this  analysis,  only  that 
portion  of  the  injected  cooling  water  which  forms  steam  and  that 
which  is  entrained  as  droplet  flow  are  taken  with  the  exhaust  gas 
flow  to  comprise  a  hypothetical  mixture  for  the  "core"  flow,  and 
only  this  flow  is  considered  pertinent  to  the  tube  bank  drag  pro¬ 
blem.  The  tube  drag  losses  associated  with  the  flow  of  the  liquid 
film  across  the  ends  of  the  spray  bars  at  the  duct  wall  are  ignored. 


Thus,  with  use  of  equation  (48)  to  express  the  pressure  loss 
coefficient,  the  pressure  drop  per  row  of  spray  bars  caused  by  tube 
drag  is  expressible  as 


AP.  = 

l 


0.020 


(l~) 


Tn-1 


(55) 


en-1 
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in  which  is  the  total  combined  exhaust  and  droplet  flux  and 

o  is  the  combined  exhaust  density  taken  equal  to  0  +  D  •  The 

s  sx  s 

total  combined  exhaust-droplet  flux  can,  by  use  of  equation  (51) , 

be  written  as 


Tn-1 


(1+Y  )  (W  +  W  .  ) 

'n-1  ex  sn-1 

A  a 
m  n-1 


(56) 


where  a  is  the  void  fraction  and  can  be  obtained,  for  example, 
from  the  empirical  correlation  of  Yagi  et  aA.  (Ref.  26) : 

0.88  0.88  0.6  -1 

(-rS  (r~)  “1°'2  +  1]  <57> 

"  e  1  1 

As  the  mixture  leaves  the  last  row  of  tubes  in  the  tube  bank,  there 
will  be  a  terminal  expansion  which  causes  an  additional  pressure 
loss  and,  again  following  the  Abramovich  analysis,  this  is  given  by 


Ap  " 


V 


eN 


(58) 


Equations  (55)  ,  (56)  ,  and  (58)  are  combined  to  give  the  pres¬ 

sure  loss  across  the  n-th  control  volume  caused  by  the  spray  bar 
drag  as 


TD 


0,020  /  S 

g  Vb 
c  o 


( 1+Y  )  (W  +  W  .  ) 

n-1  ex  sn-1 

2  2 
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(59) 


The  total  pressure  drop  across  the  n-th  control  volume  is  then  ex¬ 
pressed  as  the  sum, 

P  -  P  ,  =  -AP 1  -  TD 
n  n-1 


or 
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P  = 


P  .  -  \  I  (1+y  )  (W  +W  )V  -  (1+y  ,  )  (W  +W  .  )V  . 
n-1  g  A  '  Yn  ex  sn  n  Yn-1  ex  sn-1  n-1 
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(60) 


where  the  last  term  (the  terminal  expansion  term)  is  included  only 
if  n=N,  the  last  control  volume.  Equation  (60)  represents  the  form 
of  the  momentum  equation  to  be  used  in  the  model. 

Consider  next  an  energy  balance  written  for  the  n-th  control 
volume  which,  under  steady  flow  conditions,  (a  severe  assumption), 
gives 

Hn  .  =  (H  -  H  .  )  +  (H ,  -  H  ) 
li  sn  sn-1  dn  dn-1 


+  (H  -  H  )  +  (H  -  H  n  ) 

exn  exn-1  Ifn  lfn-1 


in  which: 
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(61) 


(62) 


(63) 

(64) 

(65) 

(66) 


The  superscript  bars  used  above  indicate  mean  values  for  the 
heat  capacities  over  the  indicated  temperature  ranges.  In  equation 
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(62)  it  has  been  assumed  that  the  liquid  droplet  evaporation  to 
form  steam  occurs  at  T  ;  this  is  consistent  with  another  assump¬ 
tion  that  must  subsequently  be  made.  Moreover,  in  equation  (65) 
the  cooling  water  temperature  rise  for  the  liquid  flowing  along  the 
duct  wall  has  been  taken  to  be  constant,  and  the  difference  in 
liquid  flow  rates  is  obtained  from  the  mass  balance  for  the  water, 
equation  (52) .  Temperature  equilibrium  between  the  droplets  and 
the  combined  exhaust  and  steam  flow  has  been  assumed  (equation  63) 
for  the  two  phase  flow  as  it  leaves  the  control  volume. 


If  the  heat  capacity  of  the  liquid  is  taken  as  a  constant  over 
the  entire  range  from  T  to  T  ,  some  simplification  of  equation 
(60)  can  be  made,  and  tKis,  when  solved  for  T  gives 
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n  L  ‘  sn-1  ps  'n-1  ex  sn-1  pi  ex  pex  )  n-1 

+  W.  .h.  (W  -  W  -)h_  -(W._  -  W._  .)C  .  (T.  T  ) 
li  li  sn  sn-1  fg  lfn  lfn-1  pi  If  cw 


+  ( v  (W  +W  )-v  ,  (W  +W 

l 'n  ex  sn  'n-1  ex  sn 

+W  C  +v(W  +W  )C  ,  I 
ex  pex  '  n'  ex  sn  pi  J 


. )  Ic  .T  /  -fw  _C 
-1  J  pi  cw  J  /  y  sn-1 


PS 


(67) 


Equations  (52) ,  (60) ,  and  (67)  comprise  the  mathematical  model 

for  the  design  problem.  The  solution  is  iterative  in  nature  and 
begins  at  the  inlet  of  the  spray  cooler,  considering  first  only  one 
row  of  spray  bars  as  the  tube  bank.  If  sufficient  reduction  of  the 
exhaust  temperature  has  not  been  achieved,  a  second  row  is  added  and 
the  problem  is  continued.  This  procedure  continues  until  the  de¬ 
sired  temperature  has  been  reached  at  the  outlet  of  the  cooler,  and 
the  pressure  drop  across  the  cooler  is  determined. 

In  this  computational  procedure,  it  is  presumed  that  the  follow¬ 
ing  information  and  system  parameters  are  either  known  or  specifiable: 


A'  V 

s, 

V 

L'  Wli' 

V  » 

Vli 

Vif,  T 

If' 

Tcw' 

Cps'  C 

Pi' 

Cpex'  hli 

Furthermore,  the  following  variables  are  assumed  to  be  known  at  the 

inlet  to  the  cooler  section:  W  ,  W  ,  P  ,  and  T  (recall  that  W 

so  ex  o  o  so 

includes  the  steam  from  both  the  ejector  and  from  the  combustion 

process,  and  that  W  is  the  "noncondensable”  exhaust  gas  flow  rate). 

GX 

Knowing  these  four  quant it ites  and  the  exhaust  gas  composition,  one 

can  obtain  n  and  V  as  follows: 
eo  o 
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in  which  y  is  the  mole  fraction  of  steam: 
so 


(68) 
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exo 

eo 
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V 
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(l-y  )p  /  R 

so  o  gex 


P  +  o 
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W  +  W 
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p  A  c/ 
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T 

o 


(69) 

(70) 


(71) 


in  which  aQ  is  unity  (equation  57)  since  the  flow  at  the  inlet  of 
the  cooler  is  assumed  to  be  only  exhaust  gas  and  steam. 


Inspection  of  equations  (52) ,  (60) ,  and  (67)  with  n=l  indicates 

the  following  quantities  are  yet  unknown:  T  ,  P..  ,  W  ,  V  ,  p  ,  and 

X  X  Si  X  6  X 

o' -j_ *  However,  if  T^ ,  P^,  and  W^are  found,  then  V^,  and  are 

readily  calculable.  Because  equation  (52)  is  used  in  equations  (60) 

and  (67) ,  one  is  left  with  two  equations  (for  T  and  P  ) ,  which  are 

n  n 

coupled  in  W  ,  and  three  unknowns  T  ,  P  ,  and  W  .  Thus  a  further 
sn  n  n  sn 

assumption  is  necessary  to  obtain  a  solution  to  the  problem,  viz . , 

that  the  relative  humidity  is  unity  for  the  combined  exhaust  flow 

as  it  leaves  the  control  volume. 


Thus,  the  following  ^wo-fold  iterative  scheme  can  be  employed: 
Assume  a  value  of  T..  =  T  and  (for  100  per  cent  relative  humidity) 

^  1  * 

obtain  the  partial  pressure  of  the  steam  P^  from  the  water  vapor 

pressure  curve  (or  calculate  it  from  a  polynomial  expression  of  the 
curve) .  The  density  is  then  calculated  as 


*  A 

p  =  P  /  R  T_ 
ysl  si  gs  1 


(72) 


Next,  assume  a  value  of  P^  =  P^  and  obtain  the  exhaust  gas  density 
from 


e  exl 


)  / 


R 

gex 


(73) 


The  total  density  is  then 
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^el  ^sl  +  ^exl 

The  steam  flow  rate  W  ,  is  calculated  as 

s  1 


W 


W  R  P  , 
ex  ex  si 


Sl  (P*  -  P*)R 

1  sl  gs 


and  the  velocity  is 


W  +  W  , 
v  =  __§x - sl_ 

1  Pel  A  “l 


(74) 


(75) 


(76) 


where  is  calculated  from  equation  (57) . 


At  this  point,  0el<  V  ,  and  have  been  obtained  based  on 

A  A 

T  and  P  .  Thus  P  can  now  be  calculated  from  equation  (60)  with 
N=l,  and  P  is  compared  with  P  .  If  P  /  P  ,  one  must  assume  a 

^  a  111  A 

different  value  for  P^  and  iterate  these  computations  until  P^  =  P^ 

to  an  acceptable  agreement.  When  P.  has  been  obtained,  equation 

(67)  is  used  to  calculate  T  ,  and  this  value  is  compared  with  T-^ . 

The  computation  is  then  iterated  on  T^  (which  necessarily  involves 

iteration  on  P^)  until  T^  is  known  to  acceptable  accuracy. 


The  calculated  value  of  T^  is  then  compared  to  the  design  temper¬ 
ature  specified  at  the  outlet  of  the  cooler  section.  If  the  exhaust 
temperature  is  too  high,  the  terminal  expansion  loss  term, 
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go  A  tv 
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(77) 


is  computed  and,  if  significant  in  magnitude,  this  term  must  be  de¬ 
leted  from  equation  (60)  ,  for  N=1  and  T ^  and  P^  re-computed  before  a 

second  row  of  spray  bars  is  added  to  the  tube  bank  and  the  problem 
solved  for  N=2 .  In  practice,  one  would  probably  not  include  this  ter¬ 
minal  expansion  loss  term  until  several  rows  of  spray  bars  had  been 
included  and  the  calculated  temperature  approaches  the  specified  de¬ 
sign  temperature  at  the  outlet  of  the  cooler  section. 
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3.3  EVALUATION  OF  THE  MODEL 

Although  development  and  refinement  of  the  model  is  continuing, 
certain  features  of  the  model  have  been  assessed  using  limited  per¬ 
formance  data  available  for  the  spray  cooler  in  the  J-5  Facility  at 
AEDC.  These  data  are  abstracted  and  given  in  Table  III.  The  rocket 
exhaust  gas  composition  given  in  Table  IV  is  taken  as  typical  for  the 
firings  which  gave  the  pressure  drop  data  given  in  Table  III.  This 
rocket  exhaust  analysis  was  used  to  obtain  the  composition  given  in 
Table  V  for  the  noncondensable  exhaust  gas,  in  this  case  the  rocket 
exhaust  gas  minus  the  steam  and  minus  the  alumina,  which  condenses 
out  of  the  gas  phase  at  about  5850  R  and  solidifies  at  about  4200  R. 
The  initial  exhaust  temperature  is  approximately  6600  R. 

The  J-5  cooler  configuration  (see  Fig.  14)  includes  eight  rows 
of  spray  bars,  only  the  first  seven  of  which  were  operative  when 
the  data  given  in  Table  III  were  obtained.  Dimensions  pertinent  to 
this  analysis  are  also  indicated  in  the  figure. 

A  momentum  equation  modified  from  equation  (60)  to  account  for 
the  unused  last  row  of  spray  bars  was  used  with  the  J-5  spray  cooler 
specifications  and  firing  data  to  test  for  consistency  between  pre¬ 
dicted  pressure  drop  and  the  measured  values.  Moreover,  some  know¬ 
ledge  of  the  behavior  of  the  droplet  entrainment  parameter  y  used  to 
specify  the  water  droplet  flow  rate  waa  obtained. 

An  iterative  procedure  was  used  to  obtain  calculated  values  of 
the  pressure  drop  versus  assumed  values  of  y  parametric  in  exhaust 
temperature  at  the  inlet  to  the  spray  cooler.  By  specifying  inlet 
temperature  and  knowing  the  inlet  pressure  and  flow  rate'  from  data 
given  in  Table  III,  one  can  evaluate  the  quantities,  such  as  exhaust 
density  and  velocity,  needed  in  the  computations.  Similarly,  knowing 
the  outlet  temperature  from  Table  III  and  using  an  estimated  pressure 
drop,  one  can  calculate  these  same  quantities  at  the  outlet  of  the 
cooler . 


At  this  point,  a  "stepped"  distribution  based  on  the  number  of 
rows  of  spray  bars  was  assumed  for  these  quantities  (viz.:  density, 
liquid  film  flow  rate,  droplet  flow  rate,  velocity,  and  the  combined 
exhaust  flow  rate)  taking  account  that  the  last  row  of  spray  bars 
was  inoperative  during  the  firings  noted  in  Table  III.  The  parameter 
y  was  also  taken  to  be  a  stepped  function  between  the  inlet  and  out¬ 
let  of  the  cooler  as 


Yn  N  VN 


n  =  1, 2 ,  .  .  .  ,N 


(78) 


in  which  n  denotes  the  n-th  control  volume  and  N  =  7  in  this  case; 
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that  is,  y  equals  y  because  the  last  row  of  spray  bars  was  not 
7  o 

used,  although  the  tube  drag  and  terminal  expansion  losses  associa¬ 
ted  with  the  last  row  were  included  in  the  computations. 

To  facilitate  these  computations,  a  polynomial  fit  was  made 
to  the  water  vapor  pressure  curve  over  a  range  of  saturation  pres¬ 
sures  from  0.75  to  7.0  psia.  This  expression, 

^  g  ^  o  ^  _ 

T*  =  0.4811P  -  7.000P  +  41.97P  +  66.29  F  (79) 

s  s  s  s 

agrees  to  within  one  percent  with  the  values  determined  experimentally 
over  the  indicated  pressure  range. 

The  calculations  were  iterated  on  the  estimated  pressure  drop 
until  the  calculated  and  estimated  values  agreed.  A  ratio  of  the 
calculated  pressure  drop  to  the  measured  pressure  drop  listed  in 
Table  III  was  then  computed.  Figures  16  through  27  show  this  ratio 
plotted  against  the  parameter  y^  parametric  in  values  of  exhaust 
inlet  temperature.  The  data  line  designation  in  the  figure  captions 
refers  to  the  corresponding  line  in  Table  III. 

In  the  computations  summarized  by  Fig.  16  through  27,  parameters 
other  than  geometric  and  thermodynamic  constants  were  taken  as: 


1)  cooling  water  supplied  to  each  row 

of  spray  bars . 10,000  gpm 

2)  average  velocity  of  spray  from  nozzles......  200  fps 

3)  average  velocity  of  liquid  film .  20  fps 


The  results  for  unit  value  of  the  pressure  drop  ratio  and  an  inlet 
temperature  of  3000  R  suggest  that  the  droplet  entrainment  parameter 
varies  between  three  and  four  based  on  these  test  firings.  The  ex¬ 
ceptional  cases  found  in  data  lines  4  and  14  (Fig.  16  and  25)  may 
reflect  measurement  errors  in  the  observed  pressure  drops  and/or  out¬ 
let  temperatures  given  in  Table  III.  The  pressure  drop  recorded  in 
data  line  14  appears  to  be  too  small  when  compared  with  similar  firing 
data  given  in  lines  15  and  16,  for  example.  In  cases  when  the  reported 

pressure  is  low,  the  y  value  corresponding  to  a  unit  value  of  the 
ratio  (calculated  pressure  drop  equal  to  experimental  pressure  drop) 
will  be  high.  This  is  necessarily  so  because  more  droplet  entrainment 
is  required  to  give  a  pressure  rise  sufficient  to  offset  more  of  the 
tube  drag  loss;  thus,  a  higher  y  value  is  obtained.  The  reported 
outlet  temperature  in  data  line  4  is  low  compared  to  similar  firings 
for  which  all  other  data  except  outlet  temperature  are  comparable  to 
line  4. 
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A  tentative  correlation  given  as  Fig.  28  is  proposed  to  esti¬ 
mate  values  of  the  droplet  entrainment  parameter  as  a  function  of 
a  correlating  parameter  which  is  effectively  the  energy  loss  asso¬ 
ciated  with  the  flow  divided  by  the  mass  flow  rate  of  the  exhaust 
stream.  This  correlation  was  made  for  unit  values  of  the  pressure 
drop  ratio  and  with  an  exhaust  inlet  temperature  of  3500  R.  Again, 
and  for  the  same  reasons  noted  above,  the  results  computed  from  data 
in  lines  4  and  14  of  Table  III  fall  away  from  the  pattern  of  the 
other  points.  These  two  points  were  thus  ignored  in  obtaining  the 
y  correlation.  The  scatter  noted  in  the  data  in  Fig.  28  could  re¬ 
sult  from  inadequacies  of  the  proposed  correlation  and/or  errors  in 
the  data  reported  in  Table  III.  While  this  correlation  does  give 
some  means  of  estimating  the  liquid  droplet  carryover  from  the  spray 
cooler  chamber,  it  does  not  obviate  the  need  for  a  much  more  funda¬ 
mental  knowledge  of  liquid  droplet  entrainment  processes  to  yield 
improved  design  of  spray  coolers. 

The  correlation  for  the  droplet  entrainment  parameter  could 
be  used  to  specify  those  values  needed  in  the  computational  proce¬ 
dure  outlined  earlier  for  the  design  model.  This  would  necessarily 
introduce  a  third  level  of  iteration  into  the  computational  scheme 
because  assumed  values  of  the  correlating  parameter  used  to  obtain 
a  v  value  would  subsequently  need  to  be  verified  with  a  computed  value 
of  the  correlating  parameter. 

To  assess  the  influence  of  certain  quantities  within  the  model 
upon  the  predicted  pressure  drop,  the  data  in  line  13  were  used  as 
a  basis  to  evaluate  parametrically  the  effects  of  changes  in  the 
cooling  water  flow  rate,  droplet  spray  discharge  velocity,  liquid 
film  velocity,  and  different  values  of  the  constant  associated  with 
the  pressure  loss  coefficient.  Results  of  these  studies  are  pre¬ 
sented  in  Fig.  29  through  32.  With  the  exception  of  the  case  for 
variation  in  the  cooling  water  flow  rate  (which  is  overwhelming  in 
all  cases) ,  there  is  a  significant  change  in  the  y  value  correspond¬ 
ing  to  a  unit  value  of  the  pressure  drop  ratio.  Of  these  parameters, 
the  liquid  film  flow  rate  is  known  with  the  least  precision. 

This  evaluation  of  the  momentum  equation  demonstrates  that  the 
approach  taken  in  the  design  model  is  consistent  with  the  experi¬ 
mental  data  available  for  spray  cooler  performance.  However,  these 
tests  have  involved  only  a  portion  of  the  total  design  model,  and 
the  further  development  of  this  work  is  being  continued  by  The 
University  of  Texas. 
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SECTION  IV 

CONCLUSIONS  AND  RECOMMENDATIONS 


It  can  be  concluded  from  the  present  investigation  that  the 
current  state  of  knowledge  of  two-phase  flow  processes  does  not 
permit  the  development  of  a  highly  sophisticated,  theoretical  flow 
model  for  spray  coolers.  At  best  one  must  resort  to  a  phenomeno¬ 
logical  approach  with  empirically  based  constants.  However,  it  has 
been  demonstrated  that  this  technique  does  permit  one  to  make  esti¬ 
mates  of  cooler  performance,  although  much  additional  work  is  needed. 

This  state  of  affairs  therefore  strongly  suggests  that,  if  fur¬ 
ther  insight  into  the  physical  mechanisms  of  cooler  flows  is  to  be 
attained,  a  comprehensive  and  systematic  experimentation  program  is 
required.  Because  spray  coolers  form  a  significant  component  of  any 
rocket  testing  facility,  it  is  obvious  that  AEDC  would  obtain  bene¬ 
fit  from  the  results  of  such  a  program,  especially  if  it  is  carried 
out  in  conjunction  with  the  normal  testing  activities  of  the  center. 
In  this  way  a  considerable  body  of  data  can  be  accumulated  in  a  re¬ 
latively  short  time  while  avoiding  any  difficulties  of  exhaust  gas 
simulation  and  modeling.  Should  such  an  on-site  spray  cooler  de¬ 
velopment  program  be  contemplated,  there  are  a  .number  of  considera¬ 
tions  which  suggest  themselves. 

Because  scaling  laws  for  two-phase  flow  are  not  well  developed, 
it  is  desirable  that  the  instrumented  test  cooler  be  as  large  as 
possible.  Short  of  running  tests  in  a  full-size  device,  it  is 
believed  that  the  test  cooler  should  be  at  least  one-fifth  to  one- 
fourth  full  size.  It  might  be  possible  to  arrange  a  duct-with- 
cooler  of  this  size  in  parallel  with  the  primary  exhaust  ducting. 

The  embryonic  state  of  measurement  techniques  with  which  to 
study  two-phase  flow  phenomena  causes  experimental  difficulties. 
However,  even  the  simplified  approach  presented  in  this  report 
could  be  refined  with  only  the  availability  of  additional  pressure 
and  temperature  measurements  and  a  more  complete  inventory  of  the 
liquid  flows.  At  present  it  appears  that  the  entrained  droplet 
flow  rate  within  the  spray  cooler  is  obtainable  only  as  the  dif¬ 
ference  between  the  cooling  water  supplied  and  that  removed  through 
the  drain  sumps;  proper  account  being  taken  for  that  water  which  is 
vaporized  to  add  steam  to  the  exhaust  flow.  Thermodynamic  properties 
and  flow  rates  of  the  exhaust  flow  into  the  cooler  section  must  also 
be  known. 

Within  the  overall  operational  objectives  of  minimizing  pres¬ 
sure  loss  and  droplet  carryover,  while  maximizing  temperature  drop. 
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some  additional  insight  could  be  obtained  by  variation  of  the  in¬ 
jected  spray  flow  rate.  This  could  be  done  during  the  normal 
course  of  operations  or  through  computer  simulation  of  the  spray 
cooler  performance  using  the  model  presented  here.  This  latter 
approach  is  being  continued  in  work  at  The  University  of  Texas. 
Because  the  geometry  of  the  spray  bar  configuration  is  important 
and,  in  many  cases,  difficult  to  treat  analytically,  the  tube  drag 
effects  for  other  than  the  simple  configurations  will  probably  be 
best  evaluated  through  experimental  investigations.  Even  cold- 
flow  pressure  loss  measurements  would  be  useful  in  this  part  of 
the  work. 

The  model  presented  here  represents  only  a  first  step  toward  the 
improved  design  of  spray  coolers  for  rocket  test  facilities.  The 
nature  of  the  problem  is  entirely  complex  and  the  need  is  for  a 
much  improved  understanding  of  two-phase  flows  from  a  fundamental 
standpoint.  With  such  knowledge  as  will  be  gained  in  the  future, 
models  such  as  that  presented  here  can  be  refined  and  extended  to 
describe  more  fully  the  performance  of  spray  cooler  systems. 
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APPENDIXES 

I.  TABLES 

II.  ILLUSTRATIONS 


AEDC-TR.  68-127 


TABLE  I 

VALUES  OF  EXPONENTS  m,  n,  AND  CONSTANTS  CL  AND  CG  FOR  VARIOUS 
LOCKHART-MARTINELLI  COMBINATIONS  OF  VISCOUS  AND  TURBULENT  FLOWS 


Flow  Types 

v-t 

t-v 

v-v 

n 

0.2 

1.0 

0.2 

1.0 

m 

0.2 

0.2 

1.0 

1.0 

CL 

0.046* 

16. 

0.046* 

16. 

CG 

0.046* 

0.046* 

16. 

16. 

(Re) liquid 

>2000 

<1000 

>2000 

<1000 

(Re)gas 

>2000 

>2000 

<1000 

<1000 

*For  smooth 

pipes 
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TABLE  II 

TEST  CONDITIONS  USED  IN  TUBE  BANK  PRESSURE  DROP  MEASUREMENTS 


Run 

G (note  1) ^ 

m 

Re 

Pq  (abs) 

lbm/sec  ft 

lbm/sec 

(note  2) 

m.Hq 

1 

32.4 

.405 

4.56xl06 

59.8 

2 

44.1 

.550 

6.20xl06 

81.0 

3 

63.3 

.788 

8.9lxl06 

111.6 

4 

79.8 

.995 

11.20xl06 

138.1 

Note: 

1.  G  based 

on  minimum  flow 

area  between  bars 

2 .  Re  based 

on  test  section 

diameter 
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TABLE  III 

PRESSURE  DROP  DATA  FROM  J-5  SPRAY  COOLER 


Data 

Line 

Firing 

Number 

Rocket 
Exhaust 
Flow 
lbm/s  ec 

Steam 

E j  ector 

Flow 

lbm/sec 

Inlet 

Pressure 

psia 

Pressure 

Drop 

psi 

Outlet 

Temperature 

o„ 

F 

4 

3402-04 

250 

101 

4.00 

o 

1 — 1 

1 

128 

5 

0407-01 

250 

89 

6.10 

-2.00 

143 

6 

0421-01 

250 

115 

5.70 

o 

CO 

l“l 

1 

141 

8 

0421-03 

237 

169 

5.10 

-1.80 

133 

9 

0421-04 

226 

100 

4.60 

I 

M 

O 

O 

12  8 

10 

0421-05 

250 

87 

6.10 

o 

CO 

I — I 

1 

143 

11 

0421-06 

250 

135 

6.60 

-1.80 

148 

12 

0447-01 

246 

326 

7.00 

1 

M 

kD 

o 

153 

13 

0447-02 

2  46 

100 

6.00 

-2.20 

140 

14 

1528-01 

30 

133 

1.48 

-0.20 

102 

15 

1528-02 

30 

150 

1.50 

-0.35 

103 

16 

1528-03 

30 

156 

1.40 

-0.30 

100 

Note:  Cooling  water  flow  was  approximately  70,000  gpm. 
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TABLE  IV 

TYPICAL  SOLID  FUEL  ROCKET  EXHAUST  COMPOSITION* 


Component 

Weiqht  Fraction 

Mole  Fraction 

H2 

0.0170 

0.2472 

H2° 

0.0487 

0.0793 

CO 

0.3586 

0.3753 

C°2 

0.0319 

0.0213 

N2 

0.1370 

0.1434 

Al2°3 

0.3749 

0.1078 

HCl 

0.0315 

0.0254 

0.9996 

0.9997 

Molecular 

weight 

=  29.31 

Specific 

gas  constant  =  52.68 

*At  temperature  « 

i  6600°R 
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TABLE  V 

"NONCONDENSABLE”  ROCKET  EXHAUST  GAS  ANALYSIS* 


Component 

Weiqht  Fraction 

Mole  Fraction 

H2 

0.0295 

0.3041 

CO 

0.6224 

0.4618 

o 

o 

ro 

0.0554 

0.0262 

N2 

0.2377 

0.1764 

HCl 

0.0547 

0.0312 

0.9997 

0.9997 

Molecular 

weight  =  20.78 

Specific 

Gas  constant  =  74.30 

*Based  on 

composition  given  in  Table 

IV  without  H^O  and  Al^O^ 
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STRATIFIELD 


WAVY 


ZZZZZ!3p* 


SLUG 


Fig.  1  Schemotic  Drowings  of  Flow  Potterns  in  Horizontol  Adiobotic 
Two-Phase  Flow 
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HOMOGENEOUS  NON¬ 
BUBBLE  HOMOGENEOUS 

OR  FROTH  BUBBLE 


SLUG 


(d)  (e)  (f)  (9)  (h) 


SEMI-ANNULAR  ANNULAR  SPRAY  MIST 

ANNULAR 


Fig.  2  Schematic  Drawings  of  Flow  Patterns  in  Adiabatic  Two-Phase 
Vertical  Upflow 
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g/X 


L  =  mass  flux  of  liquid  phase 
G  =  mass  flux  of  gas  phase 

A'-'/’gA,  )V2 

+~/ml1/3  /  v  /’l273 
v  =  surface  tension 

Fig.  3  Baker  Flow  Pattern  Chart  for  Two-Phase  Flow 
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Parameter 
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CD 

15 

3 

Eh 

4-1 

O 

C 

0 

•H 

4-1 

o 

<0 

u 

fa 


fa 


X  Parameter 


Fig.  4  Results  of  Lockhart-Mortinelli  Correlation 
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Filled  with  Liquid 
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FREE  SHEAR 
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APPROXIMATE 
PROFILE  AT 
PROBE 


Fig.  11  Relation  Between  Measured  Velocity  Profile  and  Average  Velocity 
Between  Adjacent  Tubes 
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Pressure  Loss  Coefficient, 


Fig.  13  Comparison  of  Present  Theory  and  Experimental  Data  With  Results  of  Other  Investigations  of  In-line  Tube  Banks 
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Ratio  of  Calculated  to  Experimental  Pressure  Drop 
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Fig.  17  Pressure  Drap  Ratio  Versus  Entrainment  Parameter  Parametric  in  Exhaust 
Inlet  Temperature:  Data  Line  5 
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2.0  2.5  3.0  3.5  4.0 


Droplet  Entrainment  Parameter,  y 


Fig.  18  Pressure  Drop  Ratio  Versus  Entrainment  Parameter  Parametric  in  Exhoust 
Inlet  Temperature;  Data  Line  6 
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Exhaust  Inlet  Temperature  =  6000°R 


Fig.  19  Pressure  Drop  Ratio  Versus  Entrainment  Parameter  Parametric  in  Exhaust 
Inlet  Temperature:  Data  Line  8 


Ratio  of  CalculaLted  to  Experimental  pressure  D'rop 
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Fig.  20  Pressure  Drop  Ratio  Versus  Entrainment  Parameter  Parametric  in  Exhaust 
Inlet  Temperature:  Data  Line  9 


6 1 


2.5 


4.0 


4.5 


3.0  3.5 

Droplet  Entrainment  Parameter, 


Fig.  21  Pressure  Drop  Rotia  Versus  Entrainment  Parameter  Parametric  in  Exhaust 
Inlet  Temperature:  Data  Line  10 
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Fig.  22  Pressure  Drop  Ratio  Versus  Entrainment  Parameter  Parametric  in  Exhaust 
Inlet  Temperature:  Data  Line  11 


69 


Ratio  of  Calculated  to  Experimental  Pressure  Drop 
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Fig.  23  Pressure  Drop  Ratio  Versus  Entrainment  Parameter  Parametric  in  Exhoust 
Inlet  Temperature:  Data  Line  12 
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Ratio  of  Calculated  to  Experimental  Pressure  Drop 
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Fig.  24  Pressure  Drap  Ratia  Versus  Entrainment  Parameter  Parametric  in  Exhaust 
Inlet  Temperature:  Data  Line  13 
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Ratio  of  Calculated  to  Experimental  pressure  Drop 
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Fig.  25  Pressure  Drop  Ratio  Versus  Entrainment  Parameter  Parametric  in  Exhaust 
Inlet  Temperature:  Data  Line  14 
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Ratio  of  Calculated  to  Experimental  Pressure  Drop 
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Fig.  26  Pressure  Drop  Ratio  Versus  Entrainment  Parameter  Parametric  In  Exhaust 
Inlet  Temperature:  Data  Line  15 
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Ratio  of  Calculated  to  Experimental  pressure  Drop 
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Fig.  27  Pressure  Drop  Ratio  Versus  Entrainment  Parameter  Parametric  in  Exhaust 
Inlet  Temperature:  Data  Line  16 
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Droplet  Entrainment  Parameter 
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0  2  4  6  8  10 


(  A  P/yO^)  /  W  ,  Flow  Energy  Loss  /  Exhaust  Mass  Flow  Rate 

Fig.  28  Correlotion  of  Droplet  Entrainment  Porameter  with  Ratio  of  Flow  Energy 
Loss  to  Exhoust  Outlet  Moss  Flow  Rate 
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Rati' 


0.80 


I 


Constant  =  0.032 


2  3 

Droplet  Entrainment  Paramet< 


Fig.  29  Effect  of  Pressure-Loss-Coefficient  Constont  on  Pressure 
o  Function  of  Entroinment  Porometer:  Doto  Line  13 


3500  R 


T  = 
o 

V.  .  =  200  fps 

f  l 

Vjf  =  20  fps 

4 

gpm/row  =  10 


0.0205 


2  3  4 

Droplet  Entrainment  Parameter,  X 


Fig.  30  Effect  of  Liquid  Film  Velocity  on  Pressure  Drop  Rotio  os  o  Function  of  Entroinment 
Porometer:  Doto  Line  13 


Ratio  of  Calculated  to  Experimental  pressure  Drop 


Droplet  Entrainment  Parameter,  If 


Fig.  31  Effect  af  Water  Spray  Injection  Velocity  an  Pressure  Drap  Ratio  as  a  Function  af 
Entrainment  Parameter:  Data  Line  13 
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Ratio  of  Calculated  to  Experimental  Pressure  Drop 
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Fig.  32  Effect  of  Coaling  Water  Flaw  Rate  Supplied  per  Raw  of  Spray  Bars 
on  Pressure  Drop  Ratio  as  a  Function  of  Entrainment  Parameter: 
Data  Line  13 
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